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Abstract 

Differential cross sections for the production of at least four jets have been measured in 
proton-proton collisions at = 8 TeV at the Large Hadron Collider using the ATLAS 
detector. Events are selected if the four anti-fc, R = 0.4 jets with the largest transverse 
momentum (px) within the rapidity range \y\ < 2.8 are well separated (AT?™'’ > 0.65), all 
have pj > 64 GeV, and include at least one jet with px > 100 GeV. The dataset corresponds 
to an integrated luminosity of 20.3 fb“h The cross sections, corrected for detector effects, 
are compared to leading-order and next-to-leading-order calculations as a function of the jet 
momenta, invariant masses, minimum and maximum opening angles and other kinematic 
variables. 
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1. Introduction 

The production of particle jets at hadron colliders such as the Large Hadron Collider (LHC) [1] provides a 
fertile testing ground for the theory describing strong interactions, Quantum Chromodynamics (QCD). In 
QCD, jet production is interpreted as the fragmentation of quarks and gluons produced in the scattering 
process followed by their subsequent hadronisation. At high transverse momenta (pj) the scattering 
of partons can be calculated using perturbative QCD (pQCD) and experimental jet measurements are 
directly related to the scattering of quarks and gluons. The large cross sections for such processes allow 
for differential measurements in a wide kinematic range and stringent testing of the underlying theory. 

This analysis studies events where at least four jets are produced in a hard-scatter process. These events 
are of particular interest as the corresponding Feynman diagrams require several vertices even at leading- 
order (LO) in the strong coupling constant as- The current state-of-the-art theoretical predictions for such 
processes are at next-to-leading-order in as (next-to-leading-order perturbative QCD, NLO pQCD) [2, 3], 
and they have recently been combined with parton shower (PS) simulations [4]. An alternative approach 
is taken by generators which provide a matrix element (ME) for the hardest 2^2 process while the 
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rest of the jets are provided by a PS model, which implements a resummation of the leading-logarithmic 
terms (e.g. Pythia 8 [5] and Herwig-h- [6]). It is also interesting to test multi-leg (i.e., 2 —> n) LO pQCD 
generators (e.g. Sherpa [7] or MadGraph [8]), since they may provide adequate descriptions of the data 
in specific kinematic regions and have the advantage of being less computationally expensive than NLO 
calculations. 

It is interesting to note that the previous ATLAS measurement of multi-jet production at ^[s -1 TeV [9] 
indicates that predictions may differ from dafa by ~ 30% even af NLO [10]. This work explores a variefy 
of kinemafic regimes and topological disfribufions fo fesf fhe validify of QCD calculations, including fhe 
PS approximafion and fhe necessity of higher-order ME in Monfe Carlo (MC) generafors. 

Addifionally, four-jel evenfs represenf a background to many ofher processes af hadron colliders. Hence, 
fhe predicfive power of fhe QCD calculations, in particular fheir abilify fo reproduce fhe shapes of fhe 
disfribufions sfudied in fhis analysis, is of general inferesf. While searches for new phenomena in mulfi-jef 
evenfs use dafa-driven fechniques fo esfimafe fhe confribufion from QCD evenfs, as was done for example 
in ref. [11], fhese mefhods are tested in MC simulations. The accuracy of fhe fheorefical predictions 
remains fherefore imporfanf. 

Three-jef evenfs have been measured differenfially by many experimenfs. Indeed if was observafions 
of such evenfs fhaf heralded fhe discovery of fhe gluon [12-15]. More recenfly, af fhe LHC, ATLAS 
has measured fhe three-jef cross secfion differentially [16] and CMS has used fhe rafio of three to two 
jet events to measure trs [12]. Event shape variables have also been measured, showing sensitivity to 
higher-order pQCD effects [18, 19]. Multi-jet cross sections have been measured previously at CMS [20], 
ATEAS [9], CDE [21, 22] and DO [23, 24], although with smaller datasets and/or lower energy, and 
generally focussed on different observables. 

This paper presents the differential cross sections for events with at least four jets, studied as a function of 
a variety of kinematic and topological variables which include momenta, masses and angles. Events are 
selected if the four anti-k, R - 0.4 jets with the largest transverse momentum within the rapidity range 
\y\ < 2.8 are well separated, all have pT > 64 GeV, and include at least one jet with px > 100 GeV. The 
measurements are corrected for detector effects. The variables are binned in the leading jet px and the 
total invariant mass, such that different regimes and configurations can be tested. The measurements are 
sensitive to the various mass scales in an event, the presence of forward jets, or the azimuthal configuration 
of the jets - that is, one jet recoiling against three, or two recoiling against two. 

The structure of the paper is as follows. Section 2 introduces the ATEAS detector. The observables and 
phase space of interest are defined in secfion 3. The MC simulafion samples sfudied in fhis work are 
summarised in secfion 4, while fhe fheory predicfions and fheir uncerfainfies are described in secfion 5. 
The frigger, jef calibrafion and dafa cleaning are presenfed in secfion 6. The unfolding of defecfor effecls is 
defailed in secfion 7. Secfion 8 provides fhe experimenfal uncerfainfies included in fhe final disfribufions. 
Einally, fhe resulfs are shown in secfion 9 and fhe conclusions are drawn in secfion 10. 


2. The ATLAS detector 

The ATEAS experimenf [25] is a multi-purpose particle physics defecfor wifh a forward-backward sym- 
mefric cylindrical geomefry and nearly An coverage in solid angle, wifh insfrumenfafion up fo \r]\ - A. 9} 

* ATLAS uses a right-handed Cartesian coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector. The 
Z-axis is taken along the beam pipe, and the x-axis points from the IP to the centre of the LHC ring. Cylindrical coordinates (r, (p) are used in 
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The layout of the detector is based on four superconducting magnet systems, which comprise a thin 
solenoid surrounding the inner tracking detectors (ID) and a barrel and two end-cap toroids generating 
the magnetic field for a large muon spectrometer. The calorimeters are located between the ID and the 
muon system. The lead/liquid-argon (LAr) electromagnetic (EM) calorimeter is split into two regions: 
the barrel (It/I < 1.475) and the end-cap (1.375 < |t/| < 3.2). The hadronic calorimeter is divided into 
four regions: the barrel (|7/| < 0.8) and the extended barrel (0.8 < |7/| < 1.7) made of scintillator/steel, 
the end-cap (1.5 < \r]\ < 3.2) with LAr/copper modules, and the forward calorimeter (3.1 < |?/| < 4.9) 
composed of LAr/copper and LAr/tungsten modules. 

A three-level trigger system [26] is used to select events for further analysis. The first level (LI) of 
the trigger reduces the event rate to less than 75 kHz using hardware-based trigger algorithms acting 
on a subset of detector information. The second level (L2) uses fast online algorithms, while the final 
frigger sfage, called fhe Even! Lilfer (EL), uses reconstrucfion soffware wifh algorifhms similar fo fhe 
offline versions. The lasf fwo soffware-based frigger levels, referred fo collectively as fhe High-Level 
Trigger (HLT), furlher reduce fhe even! rafe fo abouf 400 Hz. 


3. Cross-section definition 

This measuremenf uses jefs reconsfrucfed wifh fhe anfi-k, algorifhm [27] wifh four-momenfum recombin- 
afion as implemenfed in fhe LasfJef package [28]. The radius paramefer is /? = 0.4. 

Cross secfions are calculafed for evenfs wifh af leasf four jefs wifhin fhe rapidify range \y\ < 2.8. Ouf of 
fhose four jefs, fhe leading one musf have pj > 100 GeV, while fhe nexf fhree musf have pj > 64 GeV. 
In addifion, fhese four jefs musf be well separafed from one anofher by A/?™" > 0.65, where A/?™" = 

inin,jg[y4](A7?,y), and ARij - (|t/,- - -i- \(f>i - This sef of criteria is also referred fo as fhe 

i*J 

‘inclusive analysis cufs’ fo differenfiate fhem from fhe cases where addifional requiremenfs are made, for 
example on fhe invarianf mass of fhe four leading jefs. The inclusive analysis cufs are mainly mofivafed 
by fhe friggers used fo selecf evenfs, described in section 6.1. 

Cross sections are measured differenfially as a function of fhe kinemafic variables defined in fable 1 ; fhe 
lisf includes momenfum variables, mass variables and angular variables. The only jefs used in all cases 
are fhe four leading ones in pj. The observables were selecfed for fheir sensifivify fo differences befween 
differenf Monfe Carlo models of QCD processes and fheir abilify fo describe fhe dynamics of fhe evenfs. 
Lor example, fhe Hj variable is oflen used fo sef fhe scale of mulfi-jef processes. The four-jel invarianf 
mass m 4 j is represenfafive of fhe largesf energy scale in fhe even! whereas m™", fhe minimum dijef 
invarianf mass, probes fhe smallesf jef-spliffing scale. The ratio TM™"/m 4 j Iherefore provides information 
abouf fhe range of energy scales relevanf fo fhe QCD calculafion. The A0™" and A//™" variables quanfify 
fhe minimum angular separation befween any fwo jefs. The azimufhal variable A0™" distinguishes evenfs 
wifh pairs of nearby jefs (which have large A(/i™") from fhe recoil of fhree jefs againsf one (leading fo small 
A0™" values). The rapidify variable A//™" works in a similar way. The and variables are 

designed fo be sensitive fo evenfs wifh forward jefs. In order fo build firsl fhe fwo jefs wifh fhe 

largesf rapidify inferval in fhe evenf are idenfified, and fhen fhe scalar sum of fhe px of fhe remaining fwo 
jefs is calculafed. 

the transverse plane, 0 being the azimuthal angle around the beam pipe. The rapidity ij is defined by j In the pseudorapidity in terms 
of the polar angle 6 as tj = -\n tan(0/2). 


4 



Name 


Definition 


Comment 


(0 

Pj 

Hj 

m4j 




A?/™" 


A0-" 






Iplj. 


central 


Transverse momentum of the /th jet 
1 7^® 

'4 \ 2^4 ' 2 ^ 1/2 


,5N - 5-' 


2 \ 1/2 


min,;^-g[i,4] {{Ei + - (p,- + p^) ) / m4j 

i*j ^ ’ 

min;,yg[i,4] (|0,- - (pji) 
i*j 

minijg[i,4] - y-^ 

i*j 

min,;y,^:g[i,4] - (f>j\ + \(f>j - (pki) 

i*J*k 

min;j,*.g[i,4] (\yi -yj\ + \yj - yk\) 

i*j*k 

Eyf^^ = max;,^£[ 1 , 4 ] (|j/; - i/y|) 




Sorted descending in px 

Scalar sum of the px of the four jets 

Invariant mass of the four jets 

Minimum invariant mass of two jets re¬ 
lative to invariant mass of four jets 

Minimum azimuthal separation of two 
jets 

Minimum rapidity separation of two 
jets 

Minimum azimuthal separation 
between any three jets 

Minimum rapidity separation between 
any three jets 

Maximum rapidity difference between 
two jets 

If A^™“ is defined by jets a and b, this 
is the scalar sum of the px of the other 
two jets, c and d (‘central’ jets) 


Table 1; Definitions of the various kinematic variables measured. Only the four jets with the largest px are con¬ 
sidered in all cases. 


Different phase-space regions are probed by binning the variables in regions defined by a lower bound on 
and m 4 j. This allows one to distinguish between the two types of topologies characterised by A0™", or 
to track the position of the leading jet with respect to the forward-backward pair in the variables. 

Table 2 summarises all the phase-space regions considered in the analysis for each of the variables. 

The resulting differential cross-section distributions are corrected for detector effects {unfolding) and 
taken to the so-called particle-jet level, or simply ‘particle level’. In the MC simulations used in the 
unfolding procedure, particle jets are built from particles with a proper lifetime r satisfying ct > 10 mm, 
including muons and neutrinos from hadron decays. The event selection described above is applied to 
particle jets to define the phase space of the unfolded results. 

Double parton interactions have not been investigated independently, so the measurement is inclusive in 
this respect. They are expected to contribute 1% or less to the results. 
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Observable 

3 

V 

> .. 

. [GeV] 

p^^'> > ... [GeV] 

m 4 j > ... [GeV] 

Ai/-- > ... 

Pj 

Hj 

m4j 

A0-" 

Ay^- 

A0™" 

Ay^- 

j^^central 

0. 

55 


100 

100 

100 

100 

100, 400, 700, 1000 

100, 400, 700, 1000 

100, 400, 700, 1000 

100, 400, 700, 1000 

100, 250, 400, 550 

100, 250, 400, 550 

500, 1000, 1500, 2000 

1,2, 3,4 


Table 2: Summary of the analysed phase-space regions, including the and bins into which each of 

the differential cross-section measurements is split (a dash indicates when the cut is not applied on a variable). The 
A/?™" and requirements, specified in the second and third columns respectively, apply to all variables. The 
observables are defined in table 1. 


Name 

Hard scattering 

lo/nlo 

PDF 

PS/UE 

Time 

Factor 

Pythia 

Pythia 8 

LO (2 ^ 2) 

CTIO 

Pythia 8 

AU2-CT10 

0.6 

Herwig-l-+ 

Herwig-h- 

LO (2 ^ 2) 

CTEQ6L1 

Herwig-h- 

UE-EE-3-CTEQ6L1 

1.4 

MadGraph+Pythia 

MadGraph 

LO (2 ^ 4) 

CTEQ6L1 

Pythia 6 

AUET2B-CTEQ6L1 

1.1 

HEJ 

HEJ 

All 

CTIO 

— 

— 

0.9 

BlackHat/Sherpa 

BlackHat/Sherpa 

NLO (2 ^ 4) 

CTIO 

— 

— 

— 

NJet/Sherpa 

NJet/Sherpa 

NLO (2 ^ 4) 

CTIO 

— 

— 

— 


Table 3: The generators used for comparison against the data are listed, together with the parton distribution 
functions (PDFs), PS algorithms, underlying event (UE) and parameter tunes. Each MC prediction is multiplied by 
a normalisation factor (last column) as described in section 5.1, except BlackHat/Sherpa and NJet/Sherpa. (t) The 
HEJ sample is based on an approximation to all orders in a^. 


4. Monte Carlo samples 

Monte Carlo samples are used to estimate experimental systematic uncertainties, deconvolve detector 
effects, and provide predictions to be compared with the data. Leading-order Monte Carlo samples are 
used for all three purposes. A set of theoretical calculations at higher orders, described in section 5, are 
also compared to the data. The full list of generators is shown in table 3. 

The samples used in the experimental studies comprise two LO 2^2 generators, Pythia8.160 [5] and 
Herwigh-h- 2.5.2 [6], and the LO multi-leg generator MadGraphS vl. 5.12 [8]. As described in the introduc¬ 
tion, LO generators are still widely used in searches for new physics, which motivates the comparison of 
their predictions to the data. 


6 





















Both Pythia and Herwig++ employ leading-logarithmic PS models matched to LO ME calculations. Py- 
THiA uses a PS algorithm based on pj ordering, while the PS model implemented in Herwig-i-i- follows an 
angular ordering. The ME calculation provided by MadGraph contains up to four outgoing partons in the 
ME. It is matched to a PS generated with Pythia 6.427 [29] using the shower A:t-jet MEM matching [30], 
where the jet-parton matching scale is set to 20 GeV. Hadronisation effects are included via the string 
model in the case of the Pythia and MadGraph samples [29], or the cluster model [31] in events simu¬ 
lated with Herwigh-h. The parton distribution functions (PDEs) used are the NEO CTIO [32] or the EO 
distributions of CTEQ6L1 [33] as shown in table 3. 

Simulations of the underlying event, including multiple parton interactions, are included in all three EO 
samples. The parameter tunes employed are the ATEAS tunes AU2 [34] and AUET2B [35] for Pythia and 
MadGraph respectively, and the Herwig-h- tune UE-EE-3 [36]. 

The multiple pp collisions within the same and neighbouring bunch crossings (pile-up) are simulated 
as additional inelastic pp collisions using Pythia 8. Einally, the interaction of particles with the ATEAS 
detector is simulated using a GEANT4-based program [37, 38]. 


5. Theoretical predictions 

The results of the measurement are compared to NEO predictions, in addition to the EO samples de¬ 
scribed in section 4. These are calculated using BlackHat/Sherpa [2, 3] and NJet/Sherpa [39, 40], and have 
been provided by their authors. They are both fixed-order calculations with no PS and no hadronisation. 
Therefore, the results are presented at the parton-jet level, that is, using jets built from partons instead of 
hadrons. Eor the high-pT phase space covered in this analysis, non-perturbative corrections are expected 
to be small [41, 42]. BlackHat performs one-loop virtual corrections using the unitarity method and on- 
shell recursion. The remaining terms of the full NEO computation are obtained with AMEGIC-i-i- [43, 44], 
part of Sherpa. NJet makes a numerical evaluation of the one-loop virtual corrections to multi-jet produc¬ 
tion in massless QCD. The Born matrix elements are evaluated with the Comix generator [45, 46] within 
Sherpa. Sherpa also performs the phase-space integration and infra-red subtraction via the Catani-Seymour 
dipole formalism. Both the BlackHat/Sherpa and NJet/Sherpa predictions use the CTIQ PDEs. 

The results are also compared to predictions provided by HEJ [47-49]. HEJ is a fully exclusive Monte 
Carlo event generator based on a perturbative cross-section calculation which approximates the hard- 
scattering ME to all orders in the strong coupling constant as for jot multiplicities of two or greater. The 
approximation is exact in the limit of large separation in rapidity between partons. The calculation uses 
the CTl® PDEs. As in the case of the NEO predictions, no PS or hadronisation are included. 

The different predictions tested are expected to display various levels of agreement in different kinematic 
configurations. The generators which combine 2 —> 2 parton matrix elements (MEs) with parton showers 
(PSs) are in principle not expected to provide a good description of the data, particularly in regions where 
the additional jets are neither soft nor collinear. A previous measurement of multi-jet cross sections at 
7 TeV by the ATEAS Collaboration [9] found that the cross section predicted by MC models typically 
disagreed with the data by (9(40%). It also found disagreements of up to 50% in the shape of the dif¬ 
ferential cross section measured as a function of or Hj. Nevertheless, there are also examples of 
exceptional cases where these calculations perform well, which adds interest to the measurement; for 
example, the same 7 TeV ATEAS paper observed that the shape of the p^^^ distribution was described by 
Pythia within just 10%. It is also interesting to test whether PSs based on an angular ordering perform 
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better in angular variables such as A0™" or A0™" than those using momentum ordering. In contrast to 
PS predictions, multi-leg matrix element calculations matched to parton showers (ME-i-PS) were seen at 
7 TeV to significantly improve the accuracy of the cross-section calculation and the shapes of the mo¬ 
mentum observables. In the present analysis, such calculations are expected to perform better in events 
with additional high-pj jets and/or large combined invariant masses of jets. This is also the type of scen¬ 
ario where HEJ is expected to perform well, since it provides an all-order description of processes with 
more than two hard jets, and it is designed to capture the hard, wide-angle emissions which a standalone 
PS approach would miss. Variables such as Ay™“, A^™" or were included in the analysis with 

this purpose in mind. Finally, the fixed-order, four-jet NLO predictions are expected to provide a better 
estimation of the cross sections than the LO calculations. Interestingly, studies at 7 TeV found that the 
NLO cross section for four-jet events was ~ 30% higher than the data [10]. 


5.1. Normalisation 

To facilitate comparison with the data, the cross sections predicted by the LO generators as well as 
HEJ are multiplied by a scale factor. The factor is such that the integrated number of events in the 
region 500 GeV < < 1.5 TeV which satisfy the inclusive analysis cuts in section 3 is equal to the 

corresponding number in data. The full set of normalisation factors is shown in table 3. No scale factor 
is ascribed to BlackHat/Sherpa and NJet/Sherpa such that the level of agreement with data can be assessed 
in light of the theoretical uncertainties, as discussed in section section 5.2. 


5.2. Theoretical uncertainties 

Theoretical uncertainties have been computed for HEJ and the NLO predictions. The sensitivity of the 
HEJ calculation to higher-order corrections was determined by the authors of the calculation by varying 
independently the renormalisation and factorisation scales by factors of ^|2, 2, 1/ V2 and 1/2 around the 
central value of HjII. The total uncertainty is the result of taking the envelope of all the variations. The 
typical size of the uncertainty is and it is not drawn on the figures for clarity. 

The central value of the renormalisation and factorisation scales used in the NJet/Sherpa and Black¬ 
Hat/Sherpa samples is also //t/2. Scale uncertainties are evaluated for NJet/Sherpa by simultaneously 
varying both scales by factors of 1/2 and 2. PDF uncertainties are obtained by reweighting the distribu¬ 
tions for all the PDF error sets using LHAPDF [50], following the recommendations from ref. [51]. The 
additional PDF sets include variations in the value of as. The sum in quadrature of the resulting scale 
and PDF variations defines fhe NLO fheorefical uncerfainfy included in fhe resulf figures in secfion 9. 
The uncerfainfy is dominafed by fhe scale componenf due fo fhe rapid drop of fhe cross secfion wifh de¬ 
creasing values of fhe renormalisafion and facforisafion scales. As a resulf, fhe uncerfainfy is significanfly 
asymmefric. 


6. Data selection and calibration 

The dafa sample used was faken during fhe period from March fo December 2012 wifh fhe LHC operating 
at a pp cenlre-of-mass energy of = 8 TeV. The application of dafa-qualify requiremenfs resulfs in an 
infegrafed luminosity of 20.3 fb“^ 
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Figure 1; Schematic of the kinematic regions in which the four different jet triggers are used, including the total 
luminosity that each of them recorded. The term 4j45 (4j65) refers to a trigger requiring at least four jets with 
Pj > 45 GeV (65 GeV), where the pf is measured at the EF level of the triggering system. The term j280 (j360) 
refers to a trigger requiring at least one jet with pj > 280 GeV (360 GeV) at the EF level. The horizontal and 
vertical axes correspond to pj^ and respectively, both calculated at the offline level (i.e., including the full 
object calibration). 


6.1. Trigger 

The events used in this analysis are selected by a combination of four jet triggers, consisting of the three 
usual levels and defined in terms of the jets produced in the event. The hardware-based LI trigger provides 
a fast decision based on the energy measured by the calorimeter. The L2 trigger performs a simple jet 
reconstruction procedure in the geometric regions identified by the LI trigger. The final decision taken 
by the EF trigger is made using jets from the region of I 77 I < 3.2, and reconstructed from topological 
clusters [52] using the anti-k? algorithm with R = 0.4. 

The four different triggers used in this paper are shown in figure 1 . Two of the triggers select events with 
at least four jets, while the remaining two select events with at least one jet at a higher pj threshold. 
Events are split into the four non-overlapping kinematic regions shown in figure 1, requiring at least four 
well-separated jets with varying pj thresholds in order to apply the corresponding trigger. This ensures 
trigger ethciencies greater than 99% for any event passing the inclusive analysis cuts. The small residual 
loss of data due to trigger inefficiency is corrected as a function of jet px using the techniques described 
in section 7. 

As noted in figure 1, three out of the four triggers only recorded a fraction of the total dataset. The con¬ 
tributions from the events selected by those three triggers are scaled by the inverse of the corresponding 
fraction. 


6.2. Jet reconstruction and calibration 

Jets are reconstructed using the anti-kf jet algorithm [27] with four-momentum recombination and radius 
parameter R = 0.4. The inputs to the jet algorithm are locally-calibrated topological clusters of calori- 
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meter cells [52], which reconstruct the three-dimensional shower topology of each particle entering the 
calorimeter. 

ATLAS has developed several jet calibration schemes [53] with different levels of complexity and differ¬ 
ent sensitivities to systematic effects. In this analysis the local cluster weighting (LCW) calibration [52] 
method is used, which classifies fopological clusfers as eifher being of elecfromagnefic or hadronic origin. 
Based on fhis classificafion, specific energy corrections are applied, improving fhe jef energy resolufion. 
The final jef energy calibrafion, generally referred fo as fhe jef energy scale, correcfs fhe average calori- 
mefer response fo reproduce fhe energy of fhe frue parficle jef. 

The jef energy scale and resolufion have been measured in pp collision dafa using fechniques described 
in references [54-56]. The effecfs of pile-up on jef energies are accounfed for by a jef-area-based cor- 
recfion [57] prior fo fhe final calibrafion, where fhe area of fhe jef is defined in rj-cf) space. Jefs are fhen 
calibrafed fo fhe hadronic energy scale using px- and zy-dependenf calibrafion facfors based on MC sim¬ 
ulations, and fheir response is correcfed based on several observables fhaf are sensifive fo fragmenfafion 
effecfs. A residual calibration is applied fo fake info accounf differences befween dafa and MC simulation 
based on a combination of several in-sifu fechniques [54]. 


6.3. Data quality criteria 

Before applying fhe selecfion fhaf defines fhe kinemafic region of inferesf, evenfs are required fo pass 
fhe frigger, as described in secfion 6.1, and fo confain a primary verfex wifh af leasf fwo fracks. Evenfs 
which confain energy deposifs in fhe calorimefer consisfenf wifh noise, or wifh incomplefe evenf dafa, are 
rejecfed. In addifion, evenfs confaining jefs pointing fo problematic calorimefer regions, or originating 
from non-collision background, cosmic rays or defector effecfs, are vetoed. These cleaning procedures 
are emulafed in fhe MC simulation used fo correcf for experimenfal effecfs, as is discussed in defail in 
secfion 7. 

No affempf is made fo exclude jefs fhaf resulf from phofons or lepfons impacfing fhe calorimefer, nor are 
fhe confribufions from such signafures correcfed for. Evenfs confaining phofons or t lepfons are expecfed 
fo confribufe less fhan 0.1% fo fhe cross secfions under sfudy. 

Disfribufions of fwo example variables and can be seen af fhe defecfor level (i.e. prior fo 
unfolding defector effecfs) in figure 2. Differenf sefs of poinfs correspond fo fhe dafa and fhe differenf MC 
generafors, which are normalised to dafa wifh fhe scale facfors indicafed in fable 3. These are consfanf 
facfors used fo facilifafe fhe comparison wifh dafa, as described in secfion 5.1. Given fhaf fhe generafors 
have only EO or even only leading-logarifhmic accuracy, fhe observed agreemenf is reasonable. 


7. Data unfolding 

Cross secfions are measured differentially in several variables, each of which is binned in p!^^ or /n 4 j. 
Each of fhe corresponding disfribufions is individually unfolded fo deconvolve defecfor effecfs such as 
inefficiencies and resolufions. The unfolding is performed using fhe Bayesian Iferafive mefhod [58, 59], 
as implemenfed in fhe RooUnfold package [60]. The algorifhm builds an unfolding mafrix sfarfing wifh 
an initial prior probabilify disfribufion faken from MC simulafion, and improves if iferafively. The mefhod 
fakes info accounf migrafions befween bins. If also correcfs fhe resulfs for fhe presence of evenfs which 
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Figure 2: Detector-level distributions of (a) and (b) pJp* for data and for example MC predictions. The MC 
predictions have passed through detector simulation. The lower panel in each plot shows the ratios of the MC 
predictions to data. For better comparison, the predictions are multiplied by the factors indicated in the legend. 


pass the selection at reconstructed-level but not at the particle level; and for detector inefficiencies, which 
have the opposite effect. The number of iterations is optimised in order to minimise the size of the 
statistical and systematic uncertainties. A lower number of iterations results in a higher dependence on 
the MC simulation, whereas higher values give larger statistical uncertainties. For the analysis presented 
in this paper, two iterations are used. 

The data are unfolded to the particle-jet level using the Pythia MC simulation to build the unfolding 
matrix. In order to construct the matrix, events are required to pass the inclusive analysis cuts at both 
the reconstructed and particle levels. The cuts require that events have at least four jets within |i/| < 2.8, 
with >100 GeV and Pj\pj\pj^ > 64 GeV. The four leading jets must in addition be separated 
by A/?™" > 0.65. For observables requiring additional kinematic cuts, these are also applied both at 
the reconstructed and particle levels. No spatial matching is performed between reconstructed-level and 
particle-level jets. 

The correlation between the observables before and after the incorporation of experimental effects tends 
to be higher for pj-based variables, such as Hj. In the case of angular variables, such as A(/>™", the 
correlation is weakened due to cases where energy resolution effects lead to re-ordering of the jet pj. 
Nevertheless, even in the case of such angular variables the entries far from the diagonal of the correlation 
matrix are significantly smaller than the diagonal elements. The binning is derived from an optimisation 
procedure such that the purity of the bins is between 70% and 90%, and the statistical uncertainty of the 
measurement is < 10%. The purity is defined as the fractional number of events per bin which do not 
migrate to other bins after the detector simulation, calculated with respect to the number of events which 
pass the particle-level cuts. 

The possible presence of biases in the unfolded spectra due to MC mismodelling of the reconstructed-level 
spectrum is evaluated using a data-driven closure test. In this study, the MC distributions are reweighted 


11 



























to match the shape of those obtained from the data, and then unfolded using the same unfolding matrix 
as for the data. A data-driven systematic uncertainty is computed by comparing the result obtained from 
this procedure and the original reweighted particle-level MC distributions. With two iterations of the 
unfolding algorithm, this systematic uncertainty is found to be negligible. 

A second unfolding uncertainty is evaluated to account for the model dependence of the efficiency with 
which both the reconstructed- and particle-level cuts are satisfied in each MC event. The systematic 
uncertainty is derived from the differences between the efficiencies calculated with Herwig-hh- and those 
calculated using Pythia. The resulting uncertainty is found to be subdominant in most cases, with typical 
sizes of 2-10%. The uncertainty is rebinned and smoothed, such that its statistical uncertainty is smaller 
than 40%. 

The statistical uncertainties are calculated with pseudo-experiments. For each pseudo-experiment, the 
data and MC distributions are reweighted event by event following a Poisson distribution centred at one. 
Each resulting Poisson replica of the data is unfolded using the corresponding fluctuated unfolding matrix. 
The random numbers for the pseudo-experiments are generated using unique seeds, following the same 
scheme used by the inclusive jet [42], dijet [61] and three-jet [16] measurements at = 7 TeV, to allow 
for possible future combination of results with the same dataset used for this analysis. 

The integral of the unfolded distributions, corresponding to the cross section in the fiducial range determ¬ 
ined by the inclusive analysis cuts, was compared for all the variables defined in the same region of phase 
space and found to agree with each other within 0.5%. 


8. Experimental uncertainties 


Several sources of experimental uncertainty are considered in this analysis. Those arising from the un¬ 
folding procedure are described in section 7. This section presents the uncertainties which arise from the 
jet energy scale (JES), jet energy resolution (JER), jet angular resolution and integrated luminosity. The 
dominant source of uncertainty in this measurement is the JES. 

The uncertainty in the JES calibration is determined in the central detector region by exploiting the trans¬ 
verse momentum balance in Z-rjet, y-t-jet or multi-jet events, which are measured in situ. The uncertain¬ 
ties in the energy of the reference object are propagated to the jet whose energy scale is being probed. 
The uncertainty in the central region is propagated to the forward region using dijet systems balanced in 
transverse momentum. The procedure is described in detail in ref. [54] . 

The total JES uncertainty is decomposed into eighteen components, which account for the uncertainty in 
the jet energy scale calibration itself, as well as uncertainties due to the pile-up subtraction procedure, 
parton flavour differences between samples, b-]&t energy scale and punch-through. Each of these uncer¬ 
tainties is incorporated as a coherent shift of the scale of the jets in the MC simulation. The energies and 
transverse momenta of all jets with px > 20 GeV and \y\ < 2.8 are varied up and down by one standard 
deviation of each uncertainty component; these components are asymmetric, i.e. the values of the up¬ 
wards and downwards variations are different. The shifts are then propagated through the unfolding. The 
unfolded distributions corresponding to the systematically varied spectra are compared one by one to the 
nominal ones, and the difference taken as the unfolded-level uncertainty due to that JES uncertainty com¬ 
ponent. The total JES uncertainty is obtained by summing all such contributions quadratically, respecting 
the sign of the variations in the event yields; that is, positive and negative event yield variations are added 
independently. 
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Statistical uncertainties on each of the JES uncertainty components are obtained by creating Poisson rep¬ 
licas of the systematically varied spectra, obtained as explained in section 7. Such statistical uncertainties 
are used to evaluate the significance of the uncertainty for each component and for each bin of all the 
differential distributions. As in the case of the unfolding uncertainty, the unfolded-level uncertainty due 
to each JES component is then rebinned and smoothed using a Gaussian kernel regression in order to get 
statistical uncertainties smaller than 40% in all bins. The typical size of the JES uncertainty is 4-15%. 

Jets may be affected by additional energy originating from pile-up interactions. This effect is corrected 
for as part of the jet energy calibration. The distributions were binned in different ranges of the average 
number of interactions per bunch crossing in order to test the possible presence of residual effects. No 
significant deviations were observed, therefore no uncertainty associated with pile-up mismodelling was 
considered beyond the pile-up uncertainty already included in the jet calibration procedure. 

The JER has been measured in data using dijet events [62], and an uncertainty was derived from the 
differences seen befween dafa and MC prediction. In general, fhe energy resolution observed in dafa is 
somewhaf worse fhan fhaf in MC simulafions. The uncertainly on fhe observables can Iherefore be evalu- 
afed by smearing fhe energy of fhe reconslrucled jefs in fhe MC simulafion. Affer applying fhis smearing 
fo fhe jefs, an allemalive unfolding mafrix is derived and used fo unfold fhe nominal MC predicfion. Then 
fhe MC disfribufion is unfolded using bofh fhe nominal and fhe smeared mafrices, and fhe difference 
befween fhe Iwo is symmelrised and laken as fhe JER syslemafic uncertainly. The lypical size of fhis 
uncerfainfy is 1-10% of fhe cross section. 

The jel angular resolution was esfimaled in MC simulafion for fhe pseudorapidify and 0 by malching 
spatially jefs al fhe reconslrucled and parlicle level, and found fo be < 2%. This is in agreemenl wilh 
in-silu measuremenls, so no syslemafic uncertainly is assigned. 

Einally, fhe uncertainly on fhe inlegraled luminosily is ±2.8%. If is derived following fhe same melhodo- 
logy as fhaf delailed in ref. [63]. 

Two examples of fhe values of fhe lolal experimenfal syslemafic uncerfainfy are shown in figure 3 for 
Iwo represenlalive variables, namely H'l and A0™". The jel energy scale and resolufion uncertainlies 
dominale in fhe majorify of bins, being larger af fhe high and low ends of fhe //t speclrum. The unfolding 
uncerfainfy is nearly as large al low values of fhe jef momenfa, and if is Iherefore an imporlanf confribulion 
in mosl of fhe A0™" bins. 


9. Results 

The various differential cross seclions measured in evenls wilh af leasl four jefs are shown in figures 4 fo 19 
for jefs reconslrucled wilh fhe anfi-k? algorilhm wilh R = 0.4. The observables used for fhe measuremenls 
are defined in fable 1. The measuremenls are performed for a wide range of jef fransverse momenfa from 
64 GeV fo several TeV, spanning Iwo orders of magnilude in pj and over seven orders of magnilude 
in cross seclion. The measured cross seclions are correcled for all defeclor effecls using fhe unfolding 
procedure described in seclion 7. The Iheorefical prediclions described in sections 4 and 5 are compared 
fo fhe unfolded resulfs. 
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Figure 3; Total systematic uncertainty in the four-jet cross section measurement for anti-A:, R = 0.4 jets as a function 
of (a) Hf and (b) A0™". In both cases the event selection corresponds to the inclusive analysis cuts, namely 

Pj^ > 64 GeV, pJp* > 100 GeV and > 0.65. Separate bands show the jet energy scale (JES) and resolution 
(JER), and the unfolding uncertainty, as well as the combined total systematic uncertainty resulting from adding in 
quadrature all the components. The total statistical uncertainty of the unfolded data spectrum is also shown. The 
luminosity uncertainty is not shown separately but is included in the total uncertainty band. 


Summary of the results The scale factors applied to LO generators (see section 5.1) are found to 
vary between 0.6 and 1.4, as shown previously in table 3. Not all generators describe the shape of 
correctly, so these scale factors should not be seen as a measure of the level of agreement between MC 
simulation and data, which may vary as a function of the cuts in ^ and The cross section predicted 
by BlackHat/Sherpa and NJet/Sherpa is larger than that measured in data, but overall the difference is 
covered by the scale and PDF uncertainties evaluated using NJet/Sherpa, with only a few exceptions. 
BlackHat/Sherpa and NJet/Sherpa give identical results within statistical uncertainties; therefore only one 
of the two (NJet/Sherpa) is discussed in the following, for simplicity. It is nevertheless interesting to 
compare experimental results with two different implementations of the same NLO pQCD calculations 
as an additional cross-check. 

In general, an excellent description of both the shape and the normalisation of the variables is given 
by NJet/Sherpa. The small differences found are covered by theoretical and statistical uncertainties in 
almost all cases; only the tails of and hint at deviations from the measured distribution. Mad- 

Graph+Pythia describes the data very well in most regions of phase space, the most significant discrepancy 
being in the slopes of p^^ and pl^^and derived variables. HEJ also provides a good description of most 
variables; the most significant discrepancy occurs for the angular variables Ai/™" and when p^ ^ is 

small. However when is large, HEJ describes better than NJet/Sherpa, which highlights one of 

the strengths of this calculation. The 2 —> 2 ME calculations matched to parton showers provide different 
levels of agreement depending on the variable studied; the only variable whose shape is reasonably well 
described by both Pythia and Herwig-i-i- is Hj. 

The following discussion is based on the results obtained after applying the particular choice of norm¬ 
alisation of the theoretical predictions as explained at the beginning of this section. NJet/Sherpa, which 


14 




































generally gives very good agreement with the data, is only discussed for those cases where some devi¬ 
ations are present. 


Momentum variables The momentum variables comprise the pj of the four leading jets and Part 
of the importance of these variables lies in their wide use in analyses, alone or as inputs to more complex 
observables. They are also interesting in themselves: it has been shown that the ratio of the NLO to 
the LO predictions is relatively flat across the spectrum with a maximum variation of approximately 
25% [10]. Perhaps surprisingly, the PS description of was found to be better than that of in the 
7 TeV multi-jet measurement published by ATLAS [9]. 

Figures 4 to 7 show the px distributions of the leading four jets. All the LO generators show a slope with 
respect to the data in the leading jet px (figure 4). The ratios of Herwig-hh- and HEJ to data are remark¬ 
ably flat above ~ 500 GeV and ~ 300 GeV respectively. MadGraph-fPythia is within the experimental 
uncertainties above ~ 300 GeV, and it is the only one with a positive slope in the ratio to data. 

The subleading jet px (figure 5) is well described by HEJ, while the LO generators show similar trends 
to those in MadGraph-i-Pythia describes both p^^^and p^Vell, as shown in figures 6 and 7. As the 
7 TeV results suggested, Pythia gives a good description of the distribution of p^^ HEJ and Herwig-hh- 
overestimate the number of events with high NJet/Sherpa shows a similar trend at high p^\ but the 
discrepancy is mostly covered by the theoretical uncertainties. shown in figure 8, exhibits features 
similar to those in p^^^ 

In summary, Pythia and Herwig-h- tend to describe the px spectrum of the leading jets with similar levels 
of agreement, whereas Pythia is better at describing p^^\ MadGraph-i-Pythia does a reasonable job for all 
of them, while HEJ and NJet/Sherpa are very good for the leading jets and less so for p^\ This could 
perhaps be improved by matching the calculations to PSs. 


Mass variables Mass variables are widely used in physics searches, and they are also sensitive to events 
with large separations between jets, which puts the HEJ and MadGraph-i-Pythia predictions to the test, as 
they are expected to be especially accurate in this regime. 

The distribution of the total invariant mass m 4 j is studied in figure 9. Pythia and MadGraph-i-Pythia describe 
the data very well. Herwig-i-i- describes the shape of the data between 1 TeV and 3-6 TeV. HEJ is mostly 
compatible with the measurement, but the ratio to data has a bump structure in the region of approximately 
1 to 2 TeV. This feature is also shared by NJet/Sherpa, but the differences with respect to the data are 
covered by the NLO uncertainties. 

The description of different splitting scales is tested in figure 10 through the variable m™"/m 4 j. This 
distribution is well described by Pythia, whereas Herwig-h- gets worse with increasing m 4 j, consistently 
overestimating the two ends of the m™"/m 4 j spectrum. MadGraph-i-Pythia provides a very good descrip¬ 
tion, with a flat ratio for all the m 4 j cuts. The HEJ prediction shows trends similar to those of Herwig-h- at 
higher values of m 4 j. These differences are covered in all cases by fhe large associafed scale uncerfainfy. 
NJet/Sherpa overesfimafes fhe number of evenfs in fhe very firsf bin, possibly due fo fhe lack of a PS, buf 
ofherwise agrees wifh fhe dafa wifhin fhe fheorefical uncerfainfies. 

Overall, MadGraph-i-Pythia provides fhe besf descripfion of mass variables. 
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Angular variables Similarly to mass variables, angular variables are able to test the description of 
events with small- and wide-angle radiation. In addition, they can also provide information on the global 
spatial distribution of the jets. High-px. large-angle radiation should be well captured by the ME-i-PS 
description of MadGraph-i-Pythia, or the all-orders approximation of HEJ - particularly the rapidity vari¬ 
ables and Ay™". PS generators are expected to perform poorly at large angles, given that 

they only contain two hard jets, and the rest is left to the soft- and collinear-enhanced PS. The fixed-order 
NLO prediction of NJet/Sherpa should provide a very good description of these variables too, as long as 
they are far from the infrared limit. This is indeed the case, and therefore no detailed comments about its 
performance are given here. 

Figure 11 compares the distributions of A0™" for different cuts in Pythia has a small downwards 

slope with respect to the data in all the ranges. MadGraph-i-Pythia also shows a small slope. The other 
generators, both LO and NLO, reproduce the data very well. Herwig-h-, in particular, provides a very 
good description of the data. 

The A0™" spectrum is shown in figure 12. The different cuts change the spatial distribution of 
the events, such that at low most events contain two jets recoiling against two, while at high the 
events where one jet recoils against three dominate. In general, the description of the data improves as p^^^ 
increases. For Pythia, the number of events where one jet recoils against three (low A0™") is significantly 

overestimated when p^^^ is low; as p^^'^ increases, the agreement improves such that the p^^^ > 1000 GeV 
region is very well described. MadGraph-i-Pythia, Herwig-h-h and HEJ are mostly in good agreement with 
data. 

Figure 13 compares the distributions of Ay™" with data. This variable is remarkably well described by 
Pythia, showing no significant trend. MadGraph-i-Pythia mostly underestimates high Ay™" values, while 
Herwig-hh- has a tendency to underestimate the low values. HEJ overestimates the number of events with 
high Ay™" values at low p^^\ but describes the data very well at larger values of p^'\ 

For the variable Ay™", presented in figure 14, the predictions provided by Pythia and Herwig-h-i- show in 
general a positive slope with respect to the data. MadGraph-pPythia reproduces the shape of the data well, 
as does HEJ for p^^ > 400 GeV. However, for smaller values of p^^^ HEJ overestimates the number of 
events at the end of the spectrum, as was the case for Ay™". 

The variable Ay™“, shown in figure 15, is very well described by HEJ in events with p^^^ > 400 GeV. 

The ratios to data in both Pythia and Herwig-h- have upwards slopes in all p^^^ bins. MadGraph-i-Pythia 
provides mostly a good description of the data, with a tendency to underestimate the extremes of the 
distribution. Interestingly, NJet/Sherpa seems to overestimate the number of events in the tail, although it 
is a statistically limited region and the comparison with BlackHat/Sherpa is not conclusive. 

In summary: NJet/Sherpa mostly agrees with the data within the uncertainties, but its ratio to data has an 
upwards trend in the tail of Ay^"’‘. HEJ provides a very good description of all angular variables for the 

region p!^'^ > 400 GeV, as expected, but shows significant discrepancies with respect to the data in all 
the rapidity variables for lower plj?^ values. It is important to keep in mind, though, that the associated 
scale uncertainties are large. MadGraph-i-Pythia describes all the data well, apart from the tail of Ay™" 
and the extreme values of Ay^“, which it underestimates. Herwig-h-h gives good descriptions of the ^ 
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variables, but fails at describing the rapidity variables. Pythia has some problems describing both three- 
jet variables, as well as These variables highlight the need to combine four-jet ME calculations 

with parton showers. 


j^pcentrai variables The variables setting a minimum forward-backward rapidity interval and measuring 
the total pt of the central jets were dehned to test the framework of HEJ. HEJ has been designed 

to describe events with two jets significantly separated in rapidity with additional, central, high-px radi¬ 
ation. These variables are also useful to describe the spatial configuration of the events, as they represent 
the forward-backward rapidity span of the jets, and whether the leading jet is among the two central ones 
or not. The NLO predictions and MadGraph-i-Pythia are also expected to be successful in this regime, 
whereas the 2 —> 2 generators with PSs are expected to be less suitable. 

The variable is studied for values of larger than 1, 2, 3 or 4, and for different cuts in p^j!\ 

In most cases, the description of the observable worsens significantly with increasing and 

Figures 16 to 19 correspond to the results for Ap^™ > 1,2,3,4. 

The generators with 2^2 MEs have problems describing the data around the threshold values where 
the contribution from different jets changes, which results in kinks in the ratio distributions. One such 
transition occurs at the value for which the leading jet is first allowed to be central. For p^^^ > 

400 GeV, this happens at > 464 GeV, at which point there is a major jump in Pythia in the second 

ratio plot of figure 16. Pythia gives in general fhe mosf discrepanf predicfion, wifh kinks in fhe rafio fo 
dafa af fhe fransifion poinfs fhaf reach differences of 70% af high p^^\ as well as global slopes. Herwig-hh- 
describes the data very well at lower Ap™^’^ values, but as Ap™“ grows its normalisation worsens, as well 

as the shape - particularly at high Pj \ 

MadGraph-pPythia provides an excellent description of the variables, especially at low p^^\ The 

agreement deteriorates at high p^^\ but it is not very much affected by the changes between different jet 
conhgurations, providing overall a very good description of the shapes. Most distributions are well de¬ 
scribed by HEJ, especially the high region; the low region shows more shape differences, 

which gef worse af large Ap^®’‘. This is compafible wifh similar observafions made in previous ATEAS 
measuremenfs performed wifh 7 TeV dafa [64], where if was also shown fhaf fhe agreemenf was signific- 
anfly improved affer inferfacing HEJ wifh a PS generafor. NJet/Sherpa has a fendency fo overesfimafe fhe 
number of evenfs wifh very low which may be correlafed wifh fhe p^^ discrepancy discussed 

earlier. If provides a very good descripfion of fhe dafa ofherwise. 

Tables 4 fo 49 in appendix A confain fhe numerical values of fhe measured differenhal cross secfions and 
their corresponding uncertainties. The quoted values correspond to the average differential cross sections 
over the bin ranges given. 
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Figure 4: The four-jet differential cross section as a function of leading jet pj compared to different theor¬ 

etical predictions: Pythia, Herwig-hh- and MadGraph-i-Pythia (top), and HEJ, NJet/Sherpa and BlackHat/Sherpa 
(bottom). For better comparison, the predictions are multiplied by the factors indicated in the legend. In each hg- 
ure, the top panel shows the full spectra and the bottom panel the ratios of the different predictions to the data. The 
solid band represents the total experimental systematic uncertainty cenbed at one. The patterned band represents 
the NLO scale and PDF uncertainties calculated from NJet/Sherpa centred at the nominal NJet/Sherpa values. The 
scale uncertainties for HEJ (not drawn) are typically Ths ratio curves are formed by the central values with 
vertical uncertainty lines resulting from the propagation of the statistical uncertainties of the predictions and those 
of the unfolded data spectrum. 
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Figure 5: Unfolded four-jet differential cross section as a function of compared to different theoretical predic¬ 
tions. The other details are as for figure 4. 


19 













































> 

<D 

o 


10 p\TLAS 

10^ 


10^ b- 


3 

CM 

O 

■Q 


10 

ID 

1 

"O 


10'^ 


10'^ 


10'^ 

CC 

2 

Q 

1.5 


1 

o 

0) 

0.5 

1- 

0 


^ ^ , E moaxa 

][s=8 TeV, 95 pb‘^ - 20.3 fb'^ -- 


70 10' 


I u 


I 

|pythia 8 (x 0.6) 
|Herwig++ (x 1.4) 

^ |MadGraph+Pythia (x 1.1) 
J -,-p^^i>100 GeV 



Total experimental 
systematic uncertainty 


10^ 2x10^ 

Pf [GeV] 


> 

<D 

o 


lO*’ H 


WATLAS 


10= ^ 


2 ^ 10 ^ i 

10^ i 


10^ i 
10 
1 

10'^ ^ 
10'^ I 
10' 


3 

■a 

"O 


^=8 TeV, 95 pb' - 20.3 fb' -= 



I Data 
|hEJ (X 0.9) 
|BlackHat/Sherpa (x 1.0) 
|NJet/Sherpa (x 1.0) 
-.-p^^>>100 GeV 


Total experimental 
systematic uncertainty 

NLO (scale ® PDF) 

‘ uncertainty 


Figure 6; Unfolded four-jet differential cross section as a function of compared to different theoretical predic¬ 
tions. The other details are as for figure 4. 
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Figure 7: Unfolded four-jet differential cross section as a function of compared to different theoretical predic¬ 
tions. The other details are as for figure 4. 
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Figure 8: Unfolded four-jet differential cross section as a function of Hi, compared to different theoretical predic¬ 
tions. The other details are as for figure 4. 
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Figure 9; Unfolded four-jet differential cross section as a function of m 4 j, compared to different theoretical predic¬ 
tions. The other details are as for figure 4. Some points in the ratio curves for NJet/Sherpa fall outside the j/-axis 
range, and thus the NLO uncertainty is shown partially, or not shown, in these particular bins. 
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Figure 10; Unfolded four-jet differential cross section as a function of compared to different theoretical 

predictions: Pythia, Herwig-i-+ and MadGraph+Pythia (top), and HEJ, NJet/Sherpa and BlackHat/Sherpa (bot¬ 
tom). For better comparison, the predictions are multiplied by the factors indicated in the legend. In each hgure, 
the left panel shows the full spectra and the right panel the ratios of the different predictions to the data, divided 
according to the selection criterion applied to m 4 j. The solid band represents the total experimental systematic 
uncertainty centred at one. The patterned band represents the NLO scale and PDF uncertainties calculated from 
NJet/Sherpa centred at the nominal NJet/Sherpa values. The scale uncertainties for HEJ (not drawn) are typically 
The ratio curves are formed by the central values and vertical uncertainty lines resulting from the propagation 
of the statistical uncertainties of the predictions and those of the unfolded data spectrum. 
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Figure 11; Unfolded four-jet differential cross section as a function of A0™", compared to different theoretical 
predictions. The other details are as for hgure 10, but here the multiple ratio plots correspond to different selection 
criteria applied to Some points in the ratio curves for NJet/Sherpa fall outside the i/-axis range, and thus the 
NLO uncertainty is shown partially, or not shown, in these particular bins. 
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Figure 12: Unfolded four-jet differential cross section as a function of A0™", compared to different theoretical 
predictions. The other details are as for figure 11. Some points in the ratio curves for NJet/Sherpa fall outside the 
{/-axis range, and thus the NLO uncertainty is shown partially, or not shown, in these particular bins. 
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Figure 13: Unfolded four-jet differential cross section as a function of At/™", compared to different theoretical 
predictions. The other details are as for figure 11. Some points in the ratio curves for NJet/Sherpa fall outside the 
t/-axis range, and thus the NLO uncertainty is shown partially, or not shown, in these particular bins. 
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Figure 14: Unfolded four-jet differential cross section as a function of At/™", compared to different theoretical 
predictions. The other details are as for figure 11. Some points in the ratio curves for NJet/Sherpa fall outside the 
t/-axis range, and thus the NLO uncertainty is shown partially, or not shown, in these particular bins. 
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Figure 15: Unfolded four-jet differential cross section as a function of compared to different theoretical 

predictions. The other details are as for figure 11. Some points in the ratio curves for NJet/Sherpa fall outside the 
{/-axis range, and thus the NLO uncertainty is shown partially, or not shown, in these particular bins. 
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Figure 16; Unfolded four-jet differential cross section as a function of with > 1, compared to different 

theoretical predictions. The other details are as for hgure 11. 
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Figure 17; Unfolded four-jet differential cross section as a function of with > 2, compared to different 

theoretical predictions. The other details are as for hgure 11. 
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Figure 18; Unfolded four-jet differential cross section as a function of with > 3, compared to different 

theoretical predictions. The other details are as for hgure 11. 
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Figure 19; Unfolded four-jet differential cross section as a function of with > 4, compared to different 

theoretical predictions. The other details are as for figure 11. Some points in the ratio curves for NJet/Sherpa fall 
outside the ^-axis range, and thus the NLO uncertainty is shown partially, or not shown, in these particular bins. 


33 





























































































































































10. Conclusion 


This paper presents unfolded differential eross seetions of events with at least four jets in pp eollisions at 
8 TeV eentre-of-mass energy. The eross seetions are studied as a funetion of a variety of kinematie and 
topologieal variables whieh inelude momenta, masses and angles. Events are seleeted if the four anti-kf 
R = 0.4 jets with the largest transverse momentum within the rapidity range \i)\ < 2.8 are well separated 
(A/?™" > 0.65), all have px > 64 GeV, and inelude at least one jet with pi > 100 GeV. The results are 
obtained from the analysis of the full dataset eolleeted by the ATLAS deteetor at the LHC in 2012, whieh 
eorresponds to a total integrated luminosity of 20.3 fb“^ The total experimental systematie uneertainty 
is typieally of the order of 10%, and it is dominated by the jet energy seale ealibration uneertainty. 

The measurements are eompared to NLO pQCD predietions provided by BlackHat/Sherpa and NJet/Sherpa, 
as well as the all-orders ealeulation provided by HEJ. Three leading-order ealeulations are also eonsidered, 
ineluding two 2 ^ 2 PS samples (Pythia and Herwigh-h) and a multi-leg ealeulation with up to four par- 
tons in the ME matehed to a PS generated by Pythia (MadGraph-i-Pythia). 

The LO eross seetions and HEJ are normalised by fixed faetors to faeilitate the eomparison of the speetra 
in the kinematie regions of interest; these faetors vary between 0.6 and 1.4 for the different samples, 
where the MadGraph-i-Pythia and HEJ samples are the ones that need the smallest eorreetions. The NLO 
predietions, BlackHat/Sherpa and NJet/Sherpa, are almost always eompatible with the data within their 
theoretieal uneertainties, whieh are found to be large ((9(30%) at low momenta) and asymmetrie. Within 
the normalisation seheme used, MadGraph-i-Pythia also provides a good deseription of the data, as does 
HEJ, espeeially at high leading jet px- The 2 ^ 2 PS ealeulations generally deseribe the data relatively 
poorly, although they are found to provide good predietions in some partieular eases: Pythia gives a very 
good predietion of p^^^ and Ai/^™, while Herwigh-h- performs well in the azimuthal angle variables. 

Looking at the individual distributions of the differential eross seetion, the deseription of the jet momenta 
is eompatible with previous measurements of the multi-jet eross seetions. It should be noted that HEJ, 
NJet/Sherpa and BlackHat/Sherpa give a very good deseription of the distributions of the leading jets but 
show some diserepaney with the data for p^^\ Lor variables that are partieularly sensitive to wide-angle 
eonfigurations and high-px radiation, sueh as masses or angles, BlackHat/Sherpa, NJet/Sherpa and Mad- 
Graph-i-Pythia do a remarkable job overall. HEJ also provides a good deseription of the data, the main 
exeeption being that it disagrees with the rapidity measurements in events with low At high the 
predietion is very good. These measurements expose the shorteomings of 2 ^ 2 parton ME-i-PS predie¬ 
tions in a variety of seenarios and highlight the importanee of the more sophistieated ealeulations. 
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A. Tables of the measured cross sections 


Bin 

Bin edges [GeV] 

dcr/d(pV') [fb/GeV] 

Cr [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold 

^lumi 

1 

100-155 

2.62 X 10^ 

0.3 

1.4 

+9.7 

5.2 

8.4 

2.8 

2 

155-235 

1.47 X 10"^ 

0.3 

0.6 

+1.2 
-7 8 

2.8 

6.4 

2.8 

3 

235-325 

4.89 X 10^ 

0.4 

0.4 

+6.7 

6 7 

1.5 

4.3 

2.8 

4 

325-420 

1.35 X 10^ 

< 0.1 

0.3 

+6.0 

6 0 

1.2 

2.5 

2.8 

5 

420-530 

3.56 X 10^ 

0.1 

0.2 

+6.4 

6 2 

1.1 

1.6 

2.8 

6 

530-650 

9.2 X 10* 

0.2 

0.3 

+7:6 

6 9 

1.3 

1.4 

2.8 

7 

650-790 

2.26 X 10* 

0.4 

0.3 

+7.5 

1.5 

1.3 

2.8 

8 

790-950 

5.34 

0.8 

0.2 

+1.1 
-7 8 

1.7 

1.4 

2.8 

9 

950-1130 

1.19 

2.1 

0.2 

+8.6 

1.8 

1.5 

2.8 

10 

1130-1350 

2.27 X 10-* 

3.8 

0.3 

9 0 

2.0 

1.7 

2.8 

11 

1350-1630 

3.40 X 10-2 

9.1 

0.3 

+n.7 
-10 7 

2.1 

1.8 

2.8 

12 

1630^000 

6.31 X 10-“* 

21.9 

0.4 

+ 14.5 
-1T2 

2.5 

1.9 

2.8 


Table 4; Measured differential four-jet cross section for R - 0.4 jets, in bins of Pj \ along with the uncertainties in 
the measurement. The events are selected using the inclusive analysis cuts. All uncertainties are given in %. 

statistical uncertainties due to the number of data (MC simulation) events. The other columns corres¬ 
pond to the experimental systematic uncertainties arising from JES, JER, unfolding and luminosity, respectively. 


Bin 

Bin edges [GeV] 

do-/d(p''^^) [fb/GeV] 

Cr [%] 

« [%] 

mjes [%] 

mjer [%] 

^unfold [^] 

^lumi [^] 

1 

64-145 

2.85 X 10^ 

0.2 

0.8 

+K^6 

-8 1 

3.7 

6.4 

2.8 

2 

145-255 

7.43 X 10^ 

0.3 

0.4 

+7.5 
-7 1 

1.7 

4.2 

2.8 

3 

255-385 

7.61 X 10^ 

0.4 

0.3 

+6.6 

0.8 

2.5 

2.8 

4 

385-535 

9.71 X 10* 

0.2 

0.3 


0.7 

2.0 

2.8 

5 

535-715 

1.36 X 10* 

0.6 

0.3 

0.8 

2.1 

2.8 

6 

715-930 

1.88 

1.2 

0.3 

+?:l 

0.9 

2.3 

2.8 

7 

930-1175 

2.44 X 10-* 

3.1 

0.4 

+i.3 
-8 5 

1.0 

2.4 

2.8 

8 

1175-3000 

4.91 X 10-^ 

9.2 

0.4 

+ 11.4 
-10.5 

1.0 

2.4 

2.8 

Table 5 

: Measured differential four-jet cross section for R - 

0.4 jets,in 

bins of p!j\ 

along with the uncertainties in 

the measurement. The events are selected using the inclusive analysis cuts. All other details are as for table 4. 


Bin 

Bin edges [GeV] 

dcr/d(p!r7 [fb/GeV] 

6^^ [%] 

« [%] 

mjes [%] 

mjer [%] 

^unfold 

^lumi 

1 

64-120 

4.85 X 10^ 

0.2 

0.7 

+s^5 

1 5 

3.4 

5.8 

2.8 

2 

120-205 

5.87 X 10^ 

0.3 

0.6 

1.3 

2.9 

2.8 

3 

205-305 

3.0 X 10^ 

0.6 

0.7 

+r8 

0.4 

1.7 

2.8 

4 

305^25 

1.94 X 10* 

0.5 

0.8 

-6 7 

0.3 

1.9 

2.8 

5 

425-570 

1.39 

1.9 

1.1 

+7.8 
-7 3 

0.2 

2.1 

2.8 

6 

570-2000 

1.16 X 10-2 

6.4 

1.3 

+ 10.8 
-8.7 

0.2 

2.7 

2.8 


Table 6: Measured differential four-jet cross section for R - 0.4 jets, in bins of pj \ along with the uncertainties in 
the measurement. The events are selected using the inclusive analysis cuts. All other details are as for table 4. 
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Bin 

Bin edges [GeV] 

do-/d(p^^^) [fb/GeV] 

[%] 

« [%] 

mjes [%] 

mjer [%] 

^unfold [^] 

^lumi [^] 

1 

64-85 

1.18 X 10^ 

0.2 

0.7 

+8.7 

-8 1 

3.6 

7.0 

2.8 

2 

85-135 

1.45 X 10"^ 

< 0.1 

0.8 

2.9 

3.1 

2.8 

3 

135-190 

7.65 X 10^ 

0.1 

1.1 

+7.0 
-7 2 

1.4 

1.5 

2.8 

4 

190-255 

6.24 X 10' 

0.4 

2.2 

+5.8 
-5 9 

0.8 

1.9 

2.8 

5 

255-330 

5.37 

1.2 

3.4 

+5.6 
-5 8 

0.7 

2.1 

2.8 

6 

330-415 

4.58 X 10-' 

3.7 

4.3 

+6.1 
-6 8 

0.7 

2.1 

2.8 

7 

415-1500 

4.65 X 10-^ 

11.9 

6.4 

+6.5 

-7.5 

0.7 

2.1 

2.8 

Table 7: 

: Measured differential four-jet cross section for R - 

0.4 jets,in 

bins of pJp', 

along with the uncertainties in 

the measurement. The events are selected using the inclusive analysis cuts. All other details are as for table 4. 


Bin 

Bin edges [GeV] 

dcr/diMj) [fb/GeV] 

« [%] 

« [%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi [%] 

1 

290-485 

9.85 X 10" 

0.2 

1.0 

+0 

8 3 

4.0 

7.8 

2.8 

2 

485-705 

4.73 X 10^ 

0.3 

0.4 

+7:7 

-7 3 

1.7 

4.6 

2.8 

3 

705-950 

9.0 X 10^ 

0.3 

0.3 

+6.7 

-M 
-6 6 

0.7 

2.2 

2.8 

4 

950-1225 

1.7 X 10^ 

0.1 

0.2 

0.6 

1.5 

2.8 

5 

1225-1530 

3.35 X 10' 

0.2 

0.3 

+6.8 
-7 0 

0.6 

1.5 

2.8 

6 

1530-1875 

6.62 

0.6 

0.2 

+ 7.3 

0.7 

1.5 

2.8 

7 

1875-2265 

1.29 

1.1 

0.3 


0.8 

1.7 

2.8 

8 

2265-2705 

2.23 X 10-' 

2.8 

0.3 

-8 4 

0.8 

1.6 

2.8 

9 

2705-7000 

4.88 X 10-3 

6.4 

0.3 

+ 10.9 

-10.6 

1.1 

1.6 

2.8 

Table 8 

: Measured differential four-jet cross section for R = 

0.4 jets, in bins of Ht, 

along with the uncertainties in 

the measurement. The events are selected using the inclusive analysis cuts. All other details are as for table 4. 


Bin 

Bin edges [GeV] 

dcr/d(m4j) [fb/GeV] 

« [%] 

CV [%] 

mjes [%] 

mjer [%] 

^unfold 

^lumi [%] 

1 

100-545 

1.36 X 103 

0.5 

2.7 

+s:s 

7 9 

5.6 

9.5 

2.8 

2 

545-735 

4.66 X 103 

0.4 

1.4 

+§:3 
-7 9 

3.5 

6.4 

2.8 

3 

735-935 

3.55 X 103 

0.4 

1.2 

+ 8.7 

-8 3 

3.0 

5.0 

2.8 

4 

935-1150 

2.2x 103 

0.4 

1.6 

+8.7 
-9 1 

2.4 

4.4 

2.8 

5 

1150-1375 

1.22 X 103 

0.6 

1.6 

+9.8 
-9 2 

2.1 

3.9 

2.8 

6 

1375-1620 

6.09 X 10^ 

0.6 

1.7 

+ 10.0 

9 1 

2.1 

2.9 

2.8 

7 

1620-1880 

3.09 X 10^ 

0.7 

2.2 

+9.0 

7^ 

2.1 

2.6 

2.8 

8 

1880-2160 

1.37 X 10^ 

1.0 

2.2 

2.1 

2.6 

2.8 

9 

2160-2460 

5.53 X 10' 

1.4 

2.7 

-9 0 

2.1 

3.1 

2.8 

10 

2460-2780 

2.28 X 10' 

2.0 

3.0 

+8.8 
-8 7 

2.1 

4.5 

2.8 

11 

2780-3115 

7.97 

3.3 

5.0 

+9.1 
-8 8 

2.1 

5.3 

2.8 

12 

3115-3460 

2.71 

3.0 

4.9 

+ 10.2 

2.1 

5.8 

2.8 

13 

3460-3810 

9.55 X 10-' 

8.4 

5.5 

+U.6 

-9 4 

2.1 

6.0 

2.8 

14 

3810-7000 

4.66 X 10-3 

12.4 

6.9 

+ 11.6 
-9.4 

2.1 

6.0 

2.8 


Table 9: Measured differential four-jet cross section for R - 0.4 jets, in bins of m^, along with the uncertainties in 
the measurement. The events are selected using the inclusive analysis cuts. All other details are as for table 4. 
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Bin 

Bin edges 

d(r/d(m”“'/m4j) [fb] 

« [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi 

1 

0-0.04 

1.31 X 10^ 

4.7 

9.8 

+9.5 

3.3 

7.0 

2.8 

2 

0.04-0.08 

4.74 X 10^ 

0.7 

2.5 


3.3 

7.0 

2.8 

3 

0.08-0.13 

1.58 X 10^ 

0.3 

1.4 

+i.'^ 

3.3 

7.0 

2.8 

4 

0.13-0.17 

1.85 X 10^ 

0.3 

1.3 


3.3 

6.0 

2.8 

5 

0.17-0.23 

1.29 X 10^ 

0.4 

1.4 

+i:§ 

-8 1 

3.3 

4.0 

2.8 

6 

0.23-0.3 

4.69 X 10^ 

0.5 

2.5 

+7.7 
-8 9 

3.3 

4.0 

2.8 

7 

0.3-0.4 

4.21 X 10^ 

1.4 

6.2 

+ 11.0 
-4.8 

2.0 

10.2 

2.8 


Table 10; Measured differential four-jet cross section for R - 0.4 jets, in bins of OT™"/m 4 j, along with the uncer¬ 
tainties in the measurement. The events are selected using the inclusive analysis cuts, as well as m 4 j > 500 GeV. 
All other details are as for table 4. 


Bin 

Bin edges d(T/d(m™"/OT 4 j) [fb] 

«[%] 

CV [%] 

mjes [%] 

MJER [%] 

^unfold [%] 

^lumi [^] 

1 

0-0.04 

1.34 X 

10 ^ 

4.6 

9.7 

+9.7 

1.9 

10.5 

2.8 

2 

0.04-0.08 

4.01 X 

10 ^ 

0.7 

2.2 

8.9 

1.9 

10.5 

2.8 

3 

0.08-0.13 

6.79 X 

10 ^ 

0.4 

1.4 

+8.9 

104 

1.9 

5.8 

2.8 

4 

0.13-0.17 

5.35 X 

10 ^ 

0.5 

1.8 

+9.4 
-7 9 

1.9 

5.5 

2.8 

5 

0.17-0.23 

2.21 X 

10 ^ 

0.6 

1.8 

+8.2 
-9 9 

1.9 

3.1 

2.8 

6 

0.23-0.3 

3.97 X 

10 ^ 

1.0 

2.2 

+7.1 
-7 5 

1.9 

3.1 

2.8 

7 

0.3-0.4 

1.78 X 

10 ^ 

2.7 

5.9 

+6.7 

-7.5 

1.9 

3.1 

2.8 


Table 11: Measured differential four-jet cross section for R - 0.4 jets, in bins of m™''/m 4 j, along with the uncer¬ 
tainties in the measurement. The events are selected using the inclusive analysis cuts, as well as m 4 j > 1000 GeV. 
All other details are as for table 4. 


Bin 

Bin edges dcr/d(m™‘’/m 4 j) [fb] 

« [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi 

1 

0-0.04 

1.28 X 10^ 

4.6 

9.6 

+h:2 

10 3 

2.1 

4.8 

2.8 

2 

0.04-0.08 

1.67 X 10^ 

1.1 

2.2 

+0 

-11 1 

2.1 

4.8 

2.8 

3 

0.08-0.13 

1.68 X 10^ 

0.8 

2.4 

+8.3 

2.1 

4.8 

2.8 

4 

0.13-0.17 

8.81 X 10^ 

0.9 

2.5 


2.1 

4.8 

2.8 

5 

0.17-0.23 

2.63 X 10^ 

0.9 

1.9 

+^.7 
-7 7 

2.1 

5.2 

2.8 

6 

0.23-0.3 

4.04 X 10^ 

1.7 

2.6 

+9.3 
-7 1 

2.1 

5.2 

2.8 

7 

0.3-0.4 

2.03 X 10^ 

2.2 

5.9 

+8:? 

-7.1 

2.1 

3.1 

2.8 


Table 12: Measured differential four-jet cross section for R - 0.4 jets, in bins of OT™"/m 4 j, along with the uncer¬ 
tainties in the measurement. The events are selected using the inclusive analysis cuts, as well as m 4 j > 1500 GeV. 
All other details are as for table 4. 


Bin 

Bin edges 

dcr/d(m™"/'«4j) [fb] 

« [%] 

CV [%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi [^] 

1 

0-0.04 

6.15 X 10^ 

5.3 

13.0 

+ 9:9 

7 1 

1.5 

13.2 

2.8 

2 

0.04-0.08 

4.66 X 10^ 

1.6 

2.8 

+ 10.6 

104 

1.5 

13.2 

2.8 

3 

0.08-0.13 

3.55 X 10^ 

1.2 

3.3 

+ 10.6 
-104 

1.5 

9.7 

2.8 

4 

0.13-0.17 

1.29 X 10^ 

0.8 

3.0 

+8.9 
-9 0 

1.5 

4.2 

2.8 

5 

0.17-0.23 

3.64 X 10^ 

1.1 

2.5 

+8.6 
-9 0 

1.5 

4.2 

2.8 

6 

0.23-0.3 

5.61 X 10^ 

1.3 

4.0 

+8.6 

9 0 

1.5 

4.2 

2.8 

7 

0.3-0.4 

1.85 X 10^ 

5.9 

11.3 

+8.6 

-9.0 

1.5 

16.4 

2.8 


Table 13: Measured differential four-jet cross section for R - 0.4 jets, in bins of OT™"/OT 4 j, along with the uncer¬ 
tainties in the measurement. The events are selected using the inclusive analysis cuts, as well as m 4 j > 2000 GeV. 
All other details are as for table 4. 
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Bin 

Bin edges 

dcr/d(A0™") [fb] 

[%] 

« [%] 

mjes [%] 

mjer [%] 

^unfold 

^lumi 

1 

0-0.15 

2.79 X 10*’ 

0.5 

1.8 

+ 1U.6 

-6 5 

3.8 

2.1 

2.8 

2 

0.15-0.3 

2.81 X 10® 

0.5 

1.7 

+96 

7 7 

3.8 

7.7 

2.8 

3 

0.3-0.45 

2.81 X 10® 

0.5 

1.6 

+r.i 

8 6 

3.5 

6.1 

2.8 

4 

0.45-0.6 

2.87 X 10® 

0.5 

2.0 

+ 7:2 

8 5 

3.5 

6.1 

2.8 

5 

0.6-0.75 

3.12 X 10® 

0.4 

1.7 

+T.1 

3.5 

7.6 

2.8 

6 

0.75-0.9 

2.62 X 10® 

0.5 

2.3 

+i.9 

-9 4 

3.7 

7.5 

2.8 

7 

0.9-1.05 

1.95 X 10® 

0.6 

2.5 

+9.5 
-8 8 

3.7 

6.6 

2.8 

8 

1.05-1.2 

1.52 X 10® 

0.8 

3.1 

+7.8 
-8 6 

3.7 

6.6 

2.8 

9 

1.2-1.35 

9.07 X 10® 

1.0 

4.8 

+7.8 
-8 6 

3.7 

11.0 

2.8 

10 

1.35-1.6 

1.6 X 10® 

1.9 

9.1 

+7.5 

-14.4 

3.7 

11.0 

2.8 


Table 14; Measured differential four-jet cross section for R - 0.4 jets, in bins of A0™", along with the uncertainties 
in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for table 4. 


Bin 

Bin edges 

dcr/d(A0™") [fb] 

cr [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold [^] 

^lumi [^] 

1 

0-0.15 

8.98 X 

10^ 

0.2 

0.4 

6 5 

0.9 

3.3 

2.8 

2 

0.15-0.3 

8.31 X 

10^ 

0.2 

0.5 

+6.4 
-6 2 

0.9 

3.0 

2.8 

3 

0.3-0.45 

7.37 X 

10^ 

0.2 

0.5 

+6.6 
-6 2 

1.5 

2.4 

2.8 

4 

0.45-0.6 

6.61 X 

10^ 

0.3 

0.5 

+ 5:9 

-6 3 

1.5 

3.7 

2.8 

5 

0.6-0.75 

6.28 X 

10“* 

0.3 

0.6 

+6.1 

6 7 

1.3 

2.5 

2.8 

6 

0.75-0.9 

4.07 X 

10“* 

0.3 

0.6 

+6.4 

64 

1.3 

2.2 

2.8 

7 

0.9-1.05 

2.41 X 

10“* 

0.4 

0.9 

+ 5:3 

1.8 

3.6 

2.8 

8 

1.05-1.2 

1.26 X 

10“* 

0.7 

1.3 

-6 3 

1.8 

2.8 

2.8 

9 

1.2-1.35 

5.47 X 

10® 

0.9 

2.1 

+ 5:9 

5 6 

1.8 

2.8 

2.8 

10 

1.35-1.6 

1.11 X 

10® 

1.5 

3.7 

+ 6:4 

- 5.1 

1.8 

2.8 

2.8 


Table 15: Measured differential four-jet cross section for R - 0.4 jets, in bins of A0™", along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 400 GeV. All other 
details are as for table 4. 


Bin 

Bin edges 

dcr/d(A0™") [fb] 

[%] 

« [%] 

«JES [%] 

mjer [%] 

^unfold [^] 

^lumi [^] 

1 

0-0.15 

3.57 X 10® 

1.1 

0.4 

+7.4 

1.8 

2.2 

2.8 

2 

0.15-0.3 

3.37 X 10® 

1.2 

0.4 

+7.i 

-7 6 

1.3 

3.7 

2.8 

3 

0.3-0.45 

2.94 X 10® 

1.3 

0.5 

+7.1 
-7 4 

1.7 

3.3 

2.8 

4 

0.45-0.6 

2.45 X 10® 

1.5 

0.5 

+ 7.6 

1.9 

2.8 

2.8 

5 

0.6-0.75 

2.26 X 10® 

1.3 

0.6 

2.1 

2.5 

2.8 

6 

0.75-0.9 

1.37 X 10® 

2.0 

0.8 

+1.6 

-8 1 

1.6 

2.7 

2.8 

7 

0.9-1.05 

7.54 X 10® 

2.5 

0.9 

+6.8 
-9 0 

1.6 

2.7 

2.8 

8 

1.05-1.2 

3.5 X 10® 

3.4 

1.6 

+9.0 

7 5 

1.6 

5.3 

2.8 

9 

1.2-1.35 

1.4 X 10® 

5.7 

2.3 

+id.6 

64 

1.6 

3.2 

2.8 

10 

1.35-1.6 

3.72 X 10® 

9.3 

4.7 

+ 8.0 

- 6.4 

1.6 

3.2 

2.8 


Table 16; Measured differential four-jet cross section for R - 0.4 jets, in bins of A0™", along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 700 GeV. All other 
details are as for table 4. 
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Bin 

Bin edges 

dcr/d(A0™") [fb] 

« [%] 

« [%] 

«JES 

[%] 

mjer [%] 

^unfold [^] 

^lumi [^] 

1 

0-0.15 

2.6 X 10^ 

4.4 

0.4 

+8 

9 

1.9 

2.5 

2.8 

2 

0.15-0.3 

2.27 X 10^ 

4.4 

0.6 


t 

2.1 

3.1 

2.8 

3 

0.3-0.45 

2.02 X 10^ 

4.8 

0.7 



2.1 

3.6 

2.8 

4 

0.45-0.6 

1.59 X 10^ 

5.8 

0.7 

+l 

1 

Q 

2.1 

3.5 

2.8 

5 

0.6-0.75 

1.5 X 10^ 

6.2 

0.7 


6 

2.1 

3.2 

2.8 

6 

0.75-0.9 

8.7 X 10' 

6.7 

0.9 

8 

§ 

2.8 

2.5 

2.8 

7 

0.9-1.05 

5.33 X 10' 

8.8 

1.2 

+ 10.1 

8 5 

2.8 

5.4 

2.8 

8 

1.05-1.2 

2.18 X 10' 

14.4 

2.2 

+ 10.0 
-10 2 

2.8 

4.8 

2.8 

9 

1.2-1.35 

5.96 

27.4 

5.6 

+9.5 
-7 9 

4.6 

4.8 

2.8 

10 

1.35-1.6 

3.71 

29.4 

3.5 

+9.5 

-7.9 

4.6 

6.7 

2.8 


Table 17: Measured differential four-jet cross section for R - 0.4 jets, in bins of A0™", along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 1000 GeV. All other 
details are as for table 4. 


Bin 

Bin edges 

dcr/d(A0™") [fb] 

cr [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold 

^lumi 

1 

0-0.25 

1.31 X 10^ 

5.6 

8.0 

+ 9:9 

8 0 

5.2 

7.9 

2.8 

2 

0.25-0.5 

4.41 X lO"* 

3.1 

3.3 

+ 10.0 

5.2 

7.9 

2.8 

3 

0.5-0.75 

1.18 X 10^ 

2.0 

4.4 

+?(^0 

8 1 

5.2 

7.9 

2.8 

4 

0.75-1 

2.25 X 10^ 

1.4 

2.3 

+ 10.0 
-7 2 

5.2 

13.9 

2.8 

5 

1-1.25 

2.89 X 10^ 

1.3 

3.5 

+9.9 
-7 4 

4.1 

8.0 

2.8 

6 

1.25-1.5 

4.21 X 10^ 

1.2 

2.2 

+9.3 

8 0 

4.1 

7.7 

2.8 

7 

1.5-1.75 

6.48 X 10^ 

0.9 

2.4 

+§:9 
-8 6 

4.0 

6.8 

2.8 

8 

1.75-2 

9.53 X 10^ 

0.7 

2.1 

+8.4 
-7 9 

3.8 

5.6 

2.8 

9 

2-2.25 

1.62 X 10® 

0.6 

2.0 

+ 8.3 

3.7 

5.6 

2.8 

10 

2.25-2.5 

2.24 X 10® 

0.5 

1.7 


3.7 

6.6 

2.8 

11 

2.5-2.75 

3.0 X 10® 

0.4 

1.6 

+y.i 
-7 9 

3.3 

8.7 

2.8 

12 

2.75-3 

2.81 X 10® 

0.4 

1.8 

+8.6 
-8 9 

3.3 

4.1 

2.8 

13 

3-3.25 

6.3 X 10® 

0.9 

4.0 

+ 10.9 
-7.8 

3.3 

4.1 

2.8 


Table 18: Measured differential four-jet cross section for R - 0.4 jets, in bins of A0™", along with the uncertainties 
in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for table 4. 
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Bin 

Bin edges 

dcr/d(A0™") [fb] 

« [%] 

« [%] 

mjes [%] 

mjer [%] 

^unfold 

^lumi 

1 

0-0.25 

1.97 X 10^ 

1.2 

2.3 

+6.2 

6 5 

0.8 

5.7 

2.8 

2 

0.25-0.5 

6.68 X 10^ 

0.6 

1.0 

+6.2 

6 3 

0.8 

5.7 

2.8 

3 

0.5-0.75 

1.68 X 10^ 

0.4 

0.6 

+6.7 

6 2 

0.8 

3.6 

2.8 

4 

0.75-1 

2.46 X 10^ 

0.3 

0.6 

+6.2 

0.8 

3.2 

2.8 

5 

1-1.25 

2.87 X 10^ 

0.3 

0.6 

-6 2 

0.8 

4.5 

2.8 

6 

1.25-1.5 

3.07 X 10^ 

0.3 

0.5 

+6.4 
-6 2 

0.9 

4.0 

2.8 

7 

1.5-1.75 

3.13 X 10^ 

0.3 

0.6 

+6.5 

6 2 

1.1 

3.2 

2.8 

8 

1.75-2 

3.08 X 10“* 

0.3 

0.6 

+6.3 

64 

1.8 

2.5 

2.8 

9 

2-2.25 

2.79 X 10“* 

0.3 

0.6 

+6.2 

6 5 

1.6 

2.2 

2.8 

10 

2.25-2.5 

2.55 X 10“* 

0.3 

0.6 

+6.3 

1.4 

2.2 

2.8 

11 

2.5-2.75 

2.48 X 10“* 

0.4 

0.7 

-7 3 

1.6 

3.1 

2.8 

12 

2.75-3 

2.1 X 10^ 

0.4 

0.7 

+6.1 

6 2 

1.6 

3.3 

2.8 

13 

3-3.25 

5.29 X 10^ 

0.7 

1.4 

+§:5 

-7.2 

1.1 

2.0 

2.8 


Table 19; Measured differential four-jet cross section for R - 0.4 jets, in bins of A0™", along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 400 GeV. All other 
details are as for table 4. 


Bin 

Bin edges 

dcr/d(A0™") [fb] 

« [%] 

[%] 

mjes [%] 

Mjer [%] 

^unfold 

^lumi 

1 

0-0.25 

8.01 X 10' 

5.8 

2.1 

+83 

-7 1 

2.4 

5.5 

2.8 

2 

0.25-0.5 

3.1 X 10^ 

2.8 

0.9 

+r.i 

-6 7 

2.4 

5.5 

2.8 

3 

0.5-0.75 

7.61 X 10^ 

2.2 

0.7 

+7:6 

-8 7 

1.1 

2.5 

2.8 

4 

0.75-1 

1.11 X 10^ 

1.5 

0.5 

+7.0 

7 3 

1.1 

2.3 

2.8 

5 

1-1.25 

1.22 X 10^ 

1.5 

0.6 

+7.4 
-7 5 

1.2 

2.8 

2.8 

6 

1.25-1.5 

1.22 X 10^ 

1.6 

0.5 

+7.6 

7 4 

1.5 

3.2 

2.8 

7 

1.5-1.75 

1.17 X 10^ 

1.6 

0.5 

+7.7 
-7 4 

1.7 

3.6 

2.8 

8 

1.75-2 

1.07 X 10^ 

1.6 

0.6 

+8.1 

1.8 

4.4 

2.8 

9 

2-2.25 

9.19 X 10^ 

1.8 

0.6 

-8 0 

2.0 

3.3 

2.8 

10 

2.25-2.5 

8.33 X 10^ 

1.9 

0.7 

+8.5 
-8 1 

2.6 

2.4 

2.8 

11 

2.5-2.75 

8.17 X 10^ 

1.9 

0.7 

+7.6 
-7 7 

1.9 

2.2 

2.8 

12 

2.75-3 

6.6 X 10^ 

2.4 

0.7 

+7.0 
-8 0 

1.9 

1.9 

2.8 

13 

3-3.25 

1.92 X 10^ 

4.0 

1.4 

+7.1 

-6.8 

3.2 

1.9 

2.8 


Table 20; Measured differential four-jet cross section for R - 0.4 jets, in bins of A0™", along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as Pj^ > 700 GeV. All other 
details are as for table 4. 
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Bin 

Bin edges 

dcr/d(A0™") [fb] 

« [9fc] 

« m 

«JES [%] 

mjer [%] 

^unfold [^] 

^lumi [^] 

1 

0-0.25 

5.54 

21.2 

2.5 

+9D 

-7 8 

1.6 

2.4 

2.8 

2 

0.25-0.5 

2.68 X 

10 ' 

11.6 

1.3 

+7.6 

1.6 

2.4 

2.8 

3 

0.5-0.75 

5.75 X 

10 ' 

7.4 

0.9 

+§!S 

1.6 

2.3 

2.8 

4 

0.75-1 

8.98 X 

10 ' 

5.7 

0.7 


1.6 

2.3 

2.8 

5 

1-1.25 

8.83 X 

10 ' 

6.0 

0.8 

+i.7 

1.6 

3.8 

2.8 

6 

1.25-1.5 

7.41 X 

10 ' 

6.3 

0.6 

+S:9 

1.7 

3.4 

2.8 

7 

1.5-1.75 

7.68 X 

10 ' 

6.3 

0.7 

-8 4 

1.7 

3.5 

2.8 

8 

1.75-2 

6.68 X 

10 ' 

7.0 

0.7 

+9.4 
-8 3 

2.1 

3.8 

2.8 

9 

2-2.25 

5.29 X 

10 ' 

7.3 

1.0 

+ 10.1 

7 9 

3.1 

4.5 

2.8 

10 

2.25-2.5 

4.76 X 

10 ' 

7.9 

1.0 

+ 10.3 

8 I 

3.6 

3.7 

2.8 

11 

2.5-2.75 

5.53 X 

10 ' 

6.7 

1.0 

+ 10.1 

8 1 

2.8 

3.3 

2.8 

12 

2.75-3 

4.83 X 

10 ' 

8.4 

1.0 

+9.1 
-9 2 

2.8 

3.3 

2.8 

13 

3-3.25 

1.41 X 

10 ' 

15.1 

1.7 

+ 10.7 

-6.7 

4.3 

5.6 

2.8 


Table 21; Measured differential four-jet cross section for/? = 0.4 jets, in bins of A0™", along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 1000 GeV. All other 
details are as for table 4. 


Bin 

Bin edges 

dcr/dCAi/™") [fb] 


[9fc] 

mjes [%] 

mjer [%] 

^unfold 

^lumi [^] 

1 

0-0.15 

6.61 X 10^ 

0.4 

1.2 

7 2 

3.5 

7.4 

2.8 

2 

0.15-0.3 

5.03 X 10" 

0.4 

1.5 

+8.2 

3.5 

4.6 

2.8 

3 

0.3-0.45 

3.68 X 10" 

0.4 

1.5 

+i!^ 

3.5 

4.6 

2.8 

4 

0.45-0.6 

2.58 X 10" 

0.6 

2.0 

7^ 

-8 3 

3.7 

4.6 

2.8 

5 

0.6-0.75 

1.86 X 10" 

0.6 

2.6 

+ 10.8 
-8 8 

3.7 

10.2 

2.8 

6 

0.75-0.9 

9.92 X 10" 

0.9 

3.7 

+9.6 
-8 8 

3.7 

7.3 

2.8 

7 

0.9-1.05 

5.56 X 10" 

1.2 

4.4 

+9.6 
-10 2 

3.7 

7.3 

2.8 

8 

1.05-1.2 

2.08 X 10" 

1.9 

9.4 

+9.6 
-10 2 

3.7 

7.3 

2.8 

9 

1.2-1.35 

8.57 X 10^ 

3.0 

12.6 

+7.9 

-102 

6.9 

7.3 

2.8 

10 

1.35-1.5 

3.49 X 10^ 

6.3 

19.3 

+7.9 

-102 

6.9 

15.9 

2.8 

11 

1.5-2 

2.33 X 10" 

10.9 

22.6 

+7.9 

-10.2 

6.9 

15.9 

2.8 


Table 22: Measured differential four-jet cross section for R — 0.4 jets, in bins of A^™", along with the uncertainties 
in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for table 4. 


Bin Bin edges dcr/djA;/™") [fb] g [%] [%] «Je7M 


1 

0-0.15 

1.54 X 10" 

0.2 

0.3 


2 

0.15-0.3 

1.15 X 10" 

0.2 

0.3 


3 

0.3-0.45 

7.96 X 10"' 

0.2 

0.4 

7^ 

4 

0.45-0.6 

5.23 X lO"* 

0.3 

0.6 

6 3 

5 

0.6-0.75 

3.27 X lO"* 

0.4 

0.7 

+6.0 
7 0 

6 

0.75-0.9 

1.59 X lO"* 

0.5 

0.9 

+6.2 

7 

0.9-1.05 

6.68 X 10" 

0.8 

1.5 

+^ 
f, s 

8 

1.05-1.2 

2.57 X 10" 

1.4 

2.6 

+7.6 

9 

1.2-1.35 

9.3 X 10" 

2.5 

5.5 

7^ 

-2-9 

10 

1.35-1.5 

1.97 X 10" 

4.6 

8.3 

+8.3 

-Z-9 

11 

1.5-2 

8.66 

9.9 

20.2 

+8.3 

-7.0 


mjer [%] 

0.9 

0.9 

1.7 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 


^unfold [91?] 

2J 

2.5 

3.2 

3.1 

2.4 

2.8 

2.8 

2.8 

14.0 

12.9 

12.9 


^lumi [91?] 
Z8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 


Table 23: Measured differential four-jet cross section for/? = 0.4 jets, in bins of A^™", along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 400 GeV. All other 
details are as for table 4. 
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Bin 

Bin edges 

dcr/dCAi/™") [fb] 

« [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold 

^lumi [^] 

1 

0-0.15 

6.37 X 10^ 

0.8 

0.3 

+7.3 

-7 4 

1.4 

2.7 

2.8 

2 

0.15-0.3 

4.64 X 10^ 

1.1 

0.3 

+7.6 
-7 4 

1.6 

2.9 

2.8 

3 

0.3-0.45 

2.89 X 10^ 

1.4 

0.6 

+7.7 
-7 8 

1.9 

3.3 

2.8 

4 

0.45-0.6 

1.77 X 10^ 

1.8 

0.6 

+7.7 
-7 7 

1.9 

2.8 

2.8 

5 

0.6-0.75 

9.97 X 10^ 

2.4 

0.8 

8 3 

2.6 

3.3 

2.8 

6 

0.75-0.9 

3.91 X 10^ 

3.5 

1.4 

+10 1 
-9 3 

2.6 

5.8 

2.8 

7 

0.9-1.05 

1.38 X 10^ 

6.0 

2.3 

+7.5 
-9 7 

0.8 

4.6 

2.8 

8 

1.05-1.2 

5.34 X 10* 

9.5 

4.7 

+7.5 
-9 7 

0.8 

4.6 

2.8 

9 

1.2-1.35 

1.6 X 10* 

20.9 

6.6 

+7.5 
-8 8 

0.8 

4.6 

2.8 

10 

1.35-1.5 

2.82 

45.2 

17.5 

+20.3 

-8.8 

0.8 

4.6 

2.8 


Table 24; Measured differential four-jet cross section for R - 0.4 jets, in bins of A^™", along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 700 GeV. All other 
details are as for table 4. 


Bin 

Bin edges 

d(r/d(Ay™") [fb] 

KT [%] 


mjes [%] 

mjer [%] 

^unfold [^] 

^lumi [^] 

1 

0-0.15 

4.69 X 10^ 

3.2 

0.4 

+8.6 

-7 7 

1.7 

2.5 

2.8 

2 

0.15-0.3 

3.33 X 10^ 

4.1 

0.5 

+8.7 

2.5 

2.7 

2.8 

3 

0.3-0.45 

1.98 X 10^ 

5.0 

0.7 


2.1 

5.3 

2.8 

4 

0.45-0.6 

1.08 X 10^ 

6.8 

0.8 

+ B 

-9 2 

2.1 

3.9 

2.8 

5 

0.6-0.75 

4.24 X 10* 

10.3 

1.3 

+ 11.7 
-8 7 

4.1 

4.3 

2.8 

6 

0.75-0.9 

1.72 X 10* 

16.8 

2.5 

+ 13.0 
-9 5 

4.1 

3.8 

2.8 

7 

0.9-1.05 

4.62 

34.1 

4.0 

+ 12.3 
-9 6 

1.0 

3.8 

2.8 

8 

1.05-1.2 

2.17 

50.6 

7.0 

+ 12.3 
-11.3 

1.0 

3.8 

2.8 


Table 25: Measured differential four-jet cross section for/? = 0.4 jets, in bins of A^™", along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 1000 GeV. All other 
details are as for table 4. 
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Bin 

Bin edges 

dcr/dCAi/™") [fb] 

« [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold 

^lumi [^] 

1 

0-0.25 

3.61 X 10^ 

1.3 

4.6 

+1.5 

-7 1 

3.4 

23.0 

2.8 

2 

0.25-0.5 

9.65 X 10^ 

0.7 

2.6 

+7.6 
-7 3 

3.4 

6.9 

2.8 

3 

0.5-0.75 

1.64 X 10® 

0.5 

1.9 

+7.6 
-7 4 

3.4 

6.9 

2.8 

4 

0.75-1 

1.89 X 10® 

0.5 

1.9 

+8.8 
-7 9 

3.4 

6.9 

2.8 

5 

1-1.25 

1.83 X 10® 

0.5 

1.6 

+9.0 
-8 7 

2.8 

6.6 

2.8 

6 

1.25-1.5 

1.62 X 10® 

0.5 

2.3 


2.8 

6.6 

2.8 

7 

1.5-1.75 

1.42 X 10® 

0.5 

1.9 

+i.8 

2.8 

6.6 

2.8 

8 

1.75-2 

1.07 X 10® 

0.7 

2.7 

+^4 

2.8 

7.3 

2.8 

9 

2-2.25 

7.88 X 10® 

0.8 

2.8 


3.7 

7.6 

2.8 

10 

2.25-2.5 

5.6 X 10® 

1.0 

3.4 

-9 4 

6.1 

6.8 

2.8 

11 

2.5-2.75 

3.67 X 10® 

1.1 

4.3 

+9.6 
-10 7 

7.0 

6.5 

2.8 

12 

2.75-3 

2.35 X 10® 

1.4 

5.0 

+ 10.2 
-11 1 

7.0 

6.5 

2.8 

13 

3-3.25 

1.09 X 10® 

2.2 

9.9 

+ 11.0 
-11.4 

7.0 

6.6 

2.8 

14 

3.25-3.5 

7.17 X 10"^ 

2.6 

6.5 

+ 11.8 

7.0 

9.4 

2.8 

15 

3.5-3.75 

3.72 X 10^ 

3.4 

7.1 

+li:S 

7.0 

16.9 

2.8 

16 

3.75-4 

1.4 X lO'* 

5.9 

32.0 

+li:A 

-10 9 

7.0 

19.7 

2.8 

17 

4-5 

2.46 X 10® 

7.8 

23.2 

+ 12.0 

-10.9 

15.2 

19.9 

2.8 


Table 26; Measured differential four-jet cross section for R - 0.4 jets, in bins of A^™", along with the uncertainties 
in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for table 4. 


Bin 

Bin edges 

dcr/dCAi/™") [fb] 

« [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold 

^lumi [^] 

1 

0-0.25 

8.7 X 10® 

0.6 

0.8 

+6.1 

6 2 

1.7 

2.8 

2.8 

2 

0.25-0.5 

2.49 X 10^ 

0.3 

0.7 

+6.3 

6 3 

1.7 

3.0 

2.8 

3 

0.5-0.75 

4.28 X 10^ 

0.3 

0.5 

+5:9 

6 5 

1.0 

2.0 

2.8 

4 

0.75-1 

4.83 X 10^ 

0.2 

0.4 

+6.3 

5 7 

1.0 

2.4 

2.8 

5 

1-1.25 

4.34 X 10^ 

0.3 

0.4 

+ 6:2 

6 3 

1.2 

3.1 

2.8 

6 

1.25-1.5 

3.58 X 10^ 

0.3 

0.5 

+6.1 

6 3 

1.2 

3.3 

2.8 

7 

1.5-1.75 

2.69 X 10^ 

0.4 

0.6 

+6.0 

6 4 

1.4 

3.1 

2.8 

8 

1.75-2 

1.86 X 10^ 

0.4 

0.7 

+6.4 

6 9 

1.4 

3.1 

2.8 

9 

2-2.25 

1.21 X 10^ 

0.5 

0.9 

+7.1 

-71 

1.3 

4.0 

2.8 

10 

2.25-2.5 

7.08 X 10® 

0.5 

1.2 

+7.7 

1.2 

4.4 

2.8 

11 

2.5-2.75 

4.01 X 10® 

0.8 

1.7 

1 5 

1.2 

3.6 

2.8 

12 

2.75-3 

2.0 X 10® 

1.3 

2.3 


1.2 

3.2 

2.8 

13 

3-3.25 

9.06 X 10® 

1.7 

3.1 

+^.7 
~9 1 

1.2 

3.2 

2.8 

14 

3.25-3.5 

3.88 X 10® 

2.9 

4.7 

+9.Q 
-10 1 

1.2 

3.6 

2.8 

15 

3.5-3.75 

1.26 X 10® 

4.6 

9.0 

+ 11.0 

1.2 

11.1 

2.8 

16 

3.75-4 

3.24 X 10® 

10.1 

16.6 

+lii 

1.2 

30.3 

2.8 

17 

4-5 

4.08 

17.1 

37.0 

+ili 

- 11.5 

1.2 

40.7 

2.8 


Table 27: Measured differential four-jet cross section for 7? = 0.4 jets, in bins of A^™", along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 400 GeV. All other 
details are as for table 4. 
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Bin 

Bin edges 

dcr/d{Ayfp [fb] 

[%] 

« [%] 

mjes [%] 

mjer [%] 

^unfold [^] 

^lumi [^] 

1 

0-0.25 

4.06 X 10^ 

2.3 

1.1 

+6.7 

7 6 

1.5 

3.9 

2.8 

2 

0.25-0.5 

1.14 X 10^ 

1.6 

0.5 

+6.6 

-71 

1.5 

3.2 

2.8 

3 

0.5-0.75 

1.91 X 10^ 

1.4 

0.5 

+r.i 

6 9 

1.5 

3.7 

2.8 

4 

0.75-1 

2.06 X 10^ 

1.0 

0.4 

+7.5 
-7 3 

1.3 

2.3 

2.8 

5 

1-1.25 

1.71 X 10^ 

1.3 

0.5 

+7.4 
-8 0 

1.5 

2.9 

2.8 

6 

1.25-1.5 

1.28 X 10^ 

1.4 

0.6 

+7.7 
-7 8 

2.0 

2.9 

2.8 

7 

1.5-1.75 

8.57 X 10^ 

1.7 

0.7 

+^4 

2.2 

2.6 

2.8 

8 

1.75-2 

5.06 X 10^ 

2.7 

0.8 


2.4 

2.1 

2.8 

9 

2-2.25 

2.69 X 10^ 

3.4 

1.3 


2.5 

1.9 

2.8 

10 

2.25-2.5 

1.39 X 10^ 

5.0 

1.9 


2.6 

2.6 

2.8 

11 

2.5-2.75 

5.72 X 10^ 

8.0 

2.7 

-11.0 

2.6 

4.7 

2.8 

12 

2.75-3 

2.39 X 10^ 

12.0 

4.6 

+ 10.9 
-12.4 

2.6 

5.5 

2.8 

13 

3-3.25 

8.54 

16.3 

8.3 

+ 13.7 
-14.0 

2.6 

5.8 

2.8 

14 

15 

16 

3.25-3.5 

3.5-3.75 

3.75^ 

4.60 

4.97 X 10^' 

3.03 X 10' 

24.3 

80.2 

82.6 

9.2 

88.5 

81.5 

+ 14.7 

Tit:? 

iWn 

-14.1 

2.6 

2.6 

2.6 

10.9 

24.9 

30.0 

2.8 

2.8 

2.8 


Table 28: Measured differential four-jet cross section for/? = 0.4 jets, in bins of A^™", along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 700 GeV. All other 
details are as for table 4. 


Bin 

Bin edges 

dcr/d{Ayfp [fb] 

« [%] 

sz m 

mjes [%] 

MJER [%] 

^unfold [^] 

^lumi [^] 

1 

0-0.25 

3.74 X 10' 

8.2 

1.4 

+9.7 

1.5 

4.0 

2.8 

2 

0.25-0.5 

8.69 X 10' 

5.1 

0.7 

T§;^ 

-7 4 

1.5 

2.2 

2.8 

3 

0.5-0.75 

1.48 X 10^ 

4.8 

0.5 

+il 

1.5 

2.2 

2.8 

4 

0.75-1 

1.56 X 10^ 

4.4 

0.5 

Ti:? 

1.9 

2.9 

2.8 

5 

1-1.25 

1.11 X 10^ 

6.2 

0.6 

Ti:9 

8 8 

2.2 

3.8 

2.8 

6 

1.25-1.5 

8.07 X 10' 

6.5 

0.7 

+ 100 
-8 8 

2.7 

4.0 

2.8 

7 

1.5-1.75 

4.44 X 10' 

8.2 

1.0 

+ 10.8 
-8 9 

3.4 

4.2 

2.8 

8 

1.75-2 

2.48 X 10' 

10.3 

1.3 

+ 10.8 
-9 5 

3.8 

3.6 

2.8 

9 

2-2.25 

6.87 

22.3 

2.3 

+ 11.6 
-100 

3.8 

3.3 

2.8 

10 

2.25-2.5 

5.32 

23.6 

3.0 

+ 118 

-11 1 

3.6 

3.3 

2.8 

11 

2.5-2.75 

1.31 

50.7 

6.8 

+ 14.3 
-13 8 

3.5 

3.3 

2.8 

12 

2.75-3 

1.32 

39.5 

7.9 

+16.5 

-20.0 

3.4 

3.3 

2.8 


Table 29: Measured differential four-jet cross section for R - 0.4 jets, in bins of A^™", along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as ^ > 1000 GeV. All other 
details are as for table 4. 


45 













Bin 

Bin edges 

drr/d(Aj/”“) [fb] 

cr [%] 

CV [%] 

mjes [%] 

mjer [%] 

^unfold [^] 

^lumi [%] 

1 

0-0.4 

1.55 X 

10^ 

4.4 

10.2 

+S.3 

3.3 

11.8 

2.8 

2 

0.4-0.8 

1.69 X 

10^ 

1.6 

5.5 

-9 1 

3.3 

11.8 

2.8 

3 

0.8-1.2 

4.44 X 10^ 

0.9 

3.3 

+8.4 

3.3 

10.8 

2.8 

4 

1.2-1.6 

7.29 X 

10^ 

0.6 

2.6 

+M 

-7 5 

3.3 

8.4 

2.8 

5 

1.6-2 

9.77 X 

10^ 

0.5 

2.2 

+7.1 
-7 1 

3.3 

5.7 

2.8 

6 

2-2.4 

1.15 X 

10® 

0.6 

1.7 

+7.1 
-7 3 

3.4 

4.9 

2.8 

7 

2.4-2.8 

1.18X 

10® 

0.5 

1.7 

+§:5 

7 4 

3.7 

4.7 

2.8 

8 

2.8-3.2 

1.09 X 

10® 

0.6 

2.0 

43 
-8 0 

3.8 

5.1 

2.8 

9 

3.2-3.6 

9.05 X 

10® 

0.7 

2.0 

+9.5 

3.7 

5.9 

2.8 

10 

3.6-4 

6.39 X 

10® 

0.6 

2.5 

+§.l 
-9 8 

3.6 

5.4 

2.8 

11 

4-4.4 

4.55 X 

10® 

0.7 

3.3 

+ 10.4 

9 5 

3.5 

5.1 

2.8 

12 

4.4-4.8 

2.46 X 

10® 

1.0 

4.2 

+ 11.7 

10 3 

3.5 

5.1 

2.8 

13 

4.8-5.2 

9.91 X 

10^ 

1.8 

5.2 

+ 1^ 

-13 3 

3.5 

8.9 

2.8 

14 

5.2-5.6 

2.78 X 

10^ 

3.7 

10.3 

+ 13.3 
-n 3 

3.5 

8.9 

2.8 


Table 30; Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncertainties 
in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for table 4. 


Bin 

Bin edges 

dcr/d(Aj/”“) [fb] 

cr [%] 

CV [%] 

mjes [%] 

mjer [%] 

^unfold [^] 

^lumi [^] 

1 

0-0.4 

3.08 X 

10® 

6.4 

7.6 

+0.2 

1.5 

4.0 

2.8 

2 

0.4-0.8 

3.0 X 10^ 

2.2 

2.4 

5 5 

1.5 

4.0 

2.8 

3 

0.8-1.2 

8.29 X 

10^ 

1.4 

1.3 

+6.5 

1.5 

4.0 

2.8 

4 

1.2-1.6 

1.33 X 

10® 

0.9 

0.8 

5 9 

1.5 

4.0 

2.8 

5 

1.6-2 

1.73 X 

10® 

0.8 

0.7 

+6.0 
-6 0 

1.5 

4.0 

2.8 

6 

2-2.4 

1.99 X 

10® 

0.9 

0.8 

+5.9 

1.5 

4.0 

2.8 

7 

2.4-2.8 

1.92 X 

10® 

0.8 

0.7 

7^ 

-6 4 

1.5 

4.1 

2.8 

8 

2.8-3.2 

1.68 X 

10® 

0.9 

0.8 

+6.8 
-6 7 

1.5 

4.2 

2.8 

9 

3.2-3.6 

1.29 X 

10® 

1.0 

0.9 

+7.7 

1.6 

4.1 

2.8 

10 

3.6-4 

9.03 X 

10^ 

1.4 

1.1 


1.9 

4.1 

2.8 

11 

4-4.4 

4.93 X 

10^ 

1.8 

1.4 

7?:? 

-7 7 

2.0 

5.2 

2.8 

12 

4.4-4.8 

2.43 X 

10^ 

3.0 

2.3 

+7.9 
-8 3 

2.0 

6.1 

2.8 

13 

4.8-5.2 

9.11 X 

10® 

5.1 

3.0 

+9.2 

2.0 

6.4 

2.8 

14 

5.2-5.6 

1.22 X 10® 

13.8 

11.9 

7S 

-8.2 

2.0 

5.0 

2.8 


Table 31; Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 250 GeV. All other 
details are as for table 4. 
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Bin 

Bin edges 

dcr/d(Ay^p [fb] 

KT [%] 

CV [%] 

mjes [%] 

mjer [%] 

^unfold 

^lumi [^] 

1 

0-0.4 

4.07 X 10^ 

2.2 

4.2 

+53 

0.8 

3.6 

2.8 

2 

0.4-0.8 

4.4 X 10^ 

0.7 

1.0 

+6.5 

0.8 

3.6 

2.8 

3 

0.8-1.2 

1.28 X 10“* 

0.4 

0.6 

5 5 

0.8 

3.2 

2.8 

4 

1.2-1.6 

2.12 X 10“* 

0.3 

0.5 

+ 5:9 

5 1 

0.8 

2.6 

2.8 

5 

1.6-2 

2.73 X 10“* 

0.3 

0.4 

6 1 

0.8 

2.6 

2.8 

6 

2-2.4 

2.88 X 10“* 

0.2 

0.4 

6 3 

1.1 

3.2 

2.8 

7 

2.4-2.8 

2.66 X 10“* 

0.3 

0.4 

+6.3 
-6 3 

1.8 

3.1 

2.8 

8 

2.8-3.2 

2.13 X 10“* 

0.3 

0.5 

+6.4 

+l'i 

-6 9 

1.7 

2.9 

2.8 

9 

3.2-3.6 

1.46 X 10"^ 

0.4 

0.6 

1.2 

2.2 

2.8 

10 

3.6-4 

8.75 X 10^ 

0.4 

0.8 

+7.6 

1.0 

2.1 

2.8 

11 

4-4.4 

4.32 X 10^ 

0.7 

1.2 


1.0 

2.6 

2.8 

12 

4.4-4.8 

1.69 X 10^ 

1.0 

1.8 

+§:f 

-9 1 

1.0 

4.1 

2.8 

13 

4.8-5.2 

4.64 X 10^ 

2.0 

3.3 

+ 9.4 

-8 3 

1.0 

9.2 

2.8 

14 

5.2-5.6 

5.92 X 10^ 

5.4 

12.6 

+9.4 

-8.3 

1.0 

16.9 

2.8 


Table 32: Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 400 GeV. All other 
details are as for table 4. 


Bin 

Bin edges 

dcr/d(Ay^p [fb] 

« [%] 

CV [%] 

mjes [%] 

mjer [%] 

^unfold 

^lumi [^] 

1 

0-0.4 

8.16 X 10^ 

5.1 

4.5 

+5.5 

6 1 

1.9 

3.0 

2.8 

2 

0.4-0.8 

8.85 X 10^ 

1.6 

1.3 

+5.6 

1.9 

3.0 

2.8 

3 

0.8-1.2 

2.65 X 10^ 

0.8 

0.8 


1.9 

2.9 

2.8 

4 

1.2-1.6 

4.34 X 10^ 

0.7 

0.6 

-6 6 

1.7 

2.8 

2.8 

5 

1.6-2 

5.3 X 10^ 

0.6 

0.5 

+6.6 

6 7 

1.7 

2.4 

2.8 

6 

2-2.4 

5.3 X 10^ 

0.6 

0.5 

+6.9 
-6 7 

1.6 

2.6 

2.8 

7 

2.4-2.8 

4.5 X 10^ 

0.7 

0.5 

+7.2 

7 0 

1.8 

2.8 

2.8 

8 

2.8-3.2 

3.4 X 10^ 

0.7 

0.7 

+ 1.6 

2.0 

3.0 

2.8 

9 

3.2-3.6 

2.06 X 10^ 

0.9 

0.8 

-8 4 

2.0 

3.8 

2.8 

10 

3.6-4 

1.1 X 10^ 

1.2 

1.2 

+9.6 

-10 1 

1.9 

4.8 

2.8 

11 

4-4.4 

4.47 X 10^ 

2.0 

1.6 

2.1 

4.2 

2.8 

12 

4.4-4.8 

1.54 X 10^ 

3.5 

3.1 

+9.4 

-11.4 

2.4 

3.9 

2.8 

13 

4.8-5.2 

3.92 X 10^ 

6.9 

8.4 

+9.1 

-12.0 

+9.1 

-\ 1.1 

2.6 

3.9 

2.8 

14 

5.2-5.6 

5.91 

17.1 

13.9 

2.6 

3.9 

2.8 


Table 33: Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncertainties 

in the measurement. The events are selected using the inclusive analysis cuts, as well as > 550 GeV. All other 
details are as for table 4. 
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Bin 

Bin edges [GeV] ( 

dcr/d(2/7“"‘™‘) [fb/GeV] 

KT [%] 

« [%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi [^] 

1 

120-170 

1.07 X 10^ 

0.6 

2.4 

7 7 

5.6 

13.5 

2.8 

2 

170-240 

1.93 X 10^ 

0.3 

1.2 

-8 3 

4.6 

8.4 

2.8 

3 

240-315 

9.21 X 10^ 

0.4 

1.0 

+7.8 
-8 0 

1.9 

5.5 

2.8 

4 

315-395 

3.76 X 10^ 

0.4 

0.7 

+7.2 
-7 3 

0.8 

4.6 

2.8 

5 

395-480 

1.46 X 10^ 

0.5 

0.7 

+6.4 

6 5 

0.7 

4.2 

2.8 

6 

480-575 

5.75 X 10^ 

0.6 

0.5 

+6.3 

0.7 

2.6 

2.8 

7 

575-680 

2.23 X 10^ 

0.5 

0.5 

6 2 

0.7 

2.1 

2.8 

8 

680-795 

8.73 X 10' 

0.2 

0.4 

+6.5 

6 7 

0.7 

2.0 

2.8 

9 

795-930 

3.4 X 10' 

0.4 

0.5 

+6.6 
-7 0 

0.7 

2.0 

2.8 

10 

930-1085 

1.25 X 10' 

0.6 

0.5 

+6.8 
-7 2 

0.8 

2.2 

2.8 

11 

1085-1260 

4.48 

1.0 

0.7 

+7.2 

1.1 

2.5 

2.8 

12 

1260-1465 

1.56 

1.4 

0.5 


1.2 

2.6 

2.8 

13 

1465-1705 

5.07 X 10-' 

2.1 

0.5 

+^:8 
-7 9 

1.2 

2.8 

2.8 

14 

1705-1980 

1.59 X 10 ' 

3.7 

0.6 

+8.2 
-8 0 

1.3 

2.6 

2.8 

15 

1980-2300 

4.43 X 10-2 

7.0 

0.7 

+8.5 
-8 5 

1.3 

2.4 

2.8 

16 

2300-5000 

1.64 X 10-2 

11.7 

0.6 

+ 11.2 
-9.6 

1.3 

3.7 

2.8 

Table 34: Measured differential four-jet cross section for R - 0.4 jets, in bins of 

ties in the measurement. The events are selected using the inclusive analysis cuts and Ay 

are as for table 4. 

along with the uncertain- 
> 1. All other details 



Bin 

Bin edges [GeV] 

dcr/djSp^"'™') [fb/GeV] 

KT m 

[%] 

“JES [%] 

“JER [%] 

^unfold [^] 

^lumi [^] 

1 

120-170 

4.3 X 10^ 

2.5 

2.3 

+5.9 

49 

2.4 

9.2 

2.8 

2 

170-240 

1.15 X 10^ 

0.9 

0.9 

+6.4 
-5 8 

2.2 

5.8 

2.8 

3 

240-315 

1.01 X 10^ 

0.8 

0.8 

+6.9 

1.6 

4.5 

2.8 

4 

315-395 

1.48 X 10^ 

0.7 

0.7 

6 9 

1.4 

6.3 

2.8 

5 

395^80 

1.13 X 10^ 

0.6 

0.5 

+6.6 

6 6 

1.3 

6.5 

2.8 

6 

480-575 

5.73 X 10^ 

0.6 

0.5 

+6.6 

6 2 

1.0 

3.7 

2.8 

7 

575-680 

2.25 X 10^ 

0.5 

0.5 

+6.7 

0.9 

2.3 

2.8 

8 

680-795 

8.79 X 10' 

0.3 

0.4 

+^:l 

0.9 

2.0 

2.8 

9 

795-930 

3.42 X 10' 

0.4 

0.5 

-7 1 

1.1 

2.0 

2.8 

10 

930-1085 

1.26 X 10' 

0.6 

0.5 

+6.9 
-7 2 

1.4 

2.2 

2.8 

11 

1085-1260 

4.51 

1.0 

0.7 

+7.3 

1.6 

2.5 

2.8 

12 

1260-1465 

1.57 

1.4 

0.5 

-7 9 

1.5 

2.6 

2.8 

13 

1465-1705 

5.10 X 10-' 

2.1 

0.5 

+8.1 
-7 9 

1.6 

2.8 

2.8 

14 

1705-1980 

1.60 X 10 ' 

3.6 

0.6 

+8.3 

1.8 

2.6 

2.8 

15 

1980-2300 

4.45 X 10-2 

7.0 

0.7 

7i:8 

8 5 

1.5 

2.4 

2.8 

16 

2300-5000 

1.65 X 10-2 

11.7 

0.6 

+ lt3 

-9.7 

1.3 

3.7 

2.8 


Table 35; Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncer¬ 
tainties in the measurement. The events are selected using the inclusive analysis cuts and > 1, as well as 

Pj^ > 250 GeV. All other details are as for table 4. 
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Bin 

Bin edges [GeV] ( 

dcr/d(E/7“"‘™‘) [fb/GeV] 

« [%] 

[%] 

“JES [%] 

“JER [%] 

^unfold [^] 

^lumi [^] 

1 

120-170 

1.78 X 10* 

0.8 

2.5 

4-b.2 

6 6 

1.4 

3.3 

2.8 

2 

170-240 

5.29 X 10' 

0.5 

1.0 

+5.9 

1.4 

3.8 

2.8 

3 

240-315 

7.88 X 10' 

0.3 

0.7 

6 5 

1.4 

4.6 

2.8 

4 

315-395 

9.39 X 10' 

0.3 

0.7 

+6.3 

6 3 

1.4 

4.3 

2.8 

5 

395^80 

8.49 X 10' 

0.3 

0.8 

+6.4 

6 3 

1.4 

4.9 

2.8 

6 

480-575 

1.37 X 10^ 

0.2 

0.4 

+6.2 

6 3 

1.2 

4.1 

2.8 

7 

575-680 

1.09 X 10^ 

0.2 

0.4 

+6.7 

6 3 

1.1 

3.6 

2.8 

8 

680-795 

6.76 X 10' 

0.3 

0.4 

+6.5 

1.2 

3.4 

2.8 

9 

795-930 

3.39 X 10' 

0.4 

0.5 

-7 0 

1.2 

2.4 

2.8 

10 

930-1085 

1.24 X 10' 

0.6 

0.5 

+6.9 
-7 1 

1.5 

2.2 

2.8 

11 

1085-1260 

4.47 

1.0 

0.7 

+7.2 
-7 6 

1.6 

2.5 

2.8 

12 

1260-1465 

1.56 

1.4 

0.5 

+8.0 

1.5 

2.6 

2.8 

13 

1465-1705 

5.06 X 10-' 

2.1 

0.5 


1.6 

2.8 

2.8 

14 

1705-1980 

1.59 X 10-' 

3.7 

0.6 

+fl 

-7 9 

1.9 

2.6 

2.8 

15 

1980-2300 

4.41 X 10-2 

7.0 

0.7 

+8.5 
-8 4 

1.6 

2.4 

2.8 

16 

2300-5000 

1.63 X 10-2 

11.8 

0.6 

+ 11.2 

-9.6 

1.3 

3.7 

2.8 

Table 36: Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncer¬ 

tainties in the measurement. The events are selected using the inclusive analysis cuts and Aj/™’' > 1, as well as 

Pj^ > 400 GeV. All other details are as for table 4. 


Bin 

Bin edges [GeV] 

dcr/d(E/?“"«“') [fb/GeV] 

« [%] 

« [%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi 

1 

120-170 

1.94 

2.6 

3.5 

+11.7 

5 8 

10.1 

16.8 

2.8 

2 

170-240 

4.79 

1.5 

1.6 

+go 

-7 5 

5.5 

8.6 

2.8 

3 

240-315 

6.03 

1.4 

1.4 

3.2 

4.8 

2.8 

4 

315-395 

7.51 

1.1 

1.2 

+ 7:3 

-7 7 

2.2 

3.9 

2.8 

5 

395^80 

9.99 

0.9 

1.0 

+6.8 
-7 7 

1.9 

3.7 

2.8 

6 

480-575 

1.19 X 10' 

0.8 

0.9 

+ 7:1 

-7 1 

1.9 

3.7 

2.8 

7 

575-680 

1.49 X 10' 

0.6 

0.8 

+ 7.3 

1.8 

3.6 

2.8 

8 

680-795 

1.81 X 10' 

0.6 

0.5 

-1 1 

1.6 

3.3 

2.8 

9 

795-930 

1.21 X 10' 

0.6 

0.5 

+7.3 
-7 1 

1.4 

3.4 

2.8 

10 

930-1085 

8.24 

0.8 

0.6 

+7.0 

-73 

1.3 

3.3 

2.8 

11 

1085-1260 

4.41 

1.0 

0.7 

+ 7:2 

-76 

1.3 

2.9 

2.8 

12 

1260-1465 

1.54 

1.4 

0.5 

+7:§ 

-7 7 

1.4 

2.8 

2.8 

13 

1465-1705 

5.01 X 10-' 

2.1 

0.5 

+7.9 
-7 8 

1.5 

2.9 

2.8 

14 

1705-1980 

1.57 X 10-' 

3.7 

0.6 

+8.0 

1.8 

2.6 

2.8 

15 

1980-2300 

4.38 X 10-2 

7.1 

0.7 

+8.f 

-8 3 

1.5 

2.4 

2.8 

16 

2300-5000 

1.62 X 10-2 

11.8 

0.6 

+ 11.1 
-9.5 

1.2 

3.7 

2.8 


Table 37; Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncer¬ 

tainties in the measurement. The events are selected using the inclusive analysis cuts and Ai/™“ > 1, as well as 

> 550 GeV. All other details are as for table 4. 
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Bin 

Bin edges [GeV] 

dcr/d(2/?“"«“*) [fb/GeV] 

Cf [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi 

1 

120-185 

1.01 X 10^ 

0.5 

2.2 

+TU:9 

8 3 

4.6 

11.1 

2.8 

2 

185-270 

1.23 X 10^ 

0.3 

1.1 

+8.9 

3.9 

6.7 

2.8 

3 

270-365 

4.38 X 10^ 

0.4 

1.1 

+i.o 

1.7 

4.4 

2.8 

4 

365^65 

1.39 X 10^ 

0.6 

0.7 

+6.^ 

6 6 

0.8 

3.9 

2.8 

5 

465-575 

4.43 X 10^ 

0.7 

0.5 

+6.2 

6 4 

0.8 

2.5 

2.8 

6 

575-700 

1.44 X 10^ 

0.6 

0.5 

+6.5 

0.8 

2.0 

2.8 

7 

700-845 

4.56 X 10' 

0.3 

0.5 


0.8 

1.9 

2.8 

8 

845-1005 

1.46 X 10' 

0.6 

0.5 

-7 1 

0.8 

2.2 

2.8 

9 

1005-1195 

4.58 

0.9 

0.8 

+7.7 
-7 6 

0.8 

2.9 

2.8 

10 

1195-1410 

1.38 

1.6 

0.7 

+8.1 
-7 8 

0.8 

3.6 

2.8 

11 

1410-1665 

4.12 X 10-' 

2.4 

0.6 

+8.7 

0.9 

3.6 

2.8 

12 

1665-1960 

1.10 X 10-' 

4.4 

0.7 

+^:8 

1.3 

3.4 

2.8 

13 

1960-5000 

4.31 X 10-3 

7.5 

0.7 

+i.9 

+SJ 

1.6 

3.0 

2.8 

Table 38: Measured differential four-jet cross section for R - 0.4 jets, in bins of 

along with the uncertain- 


ties in the measurement. The events are selected using the inclusive analysis cuts and Ay 

> 2. All other details 


are as for table 4. 








Bin 

Bin edges [GeV] 

dcr/d(2/7“"“^) [fb/GeV] 

« [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi [^] 

1 

120-185 

3.96 X 10^ 

2.1 

1.9 

+5.7 

1.8 

8.4 

2.8 

2 

185-270 

8.47 X 10^ 

0.9 

0.8 

+y 

6 6 

1.7 

6.1 

2.8 

3 

270-365 

8.57 X 10^ 

0.9 

1.0 

+7.1 

6 8 

1.3 

7.3 

2.8 

4 

365-465 

9.37 X 10^ 

0.7 

0.5 

+7.0 

1.1 

6.9 

2.8 

5 

465-575 

4.35 X 10^ 

0.7 

0.5 

6 5 

1.1 

3.8 

2.8 

6 

575-700 

1.45 X 10^ 

0.6 

0.5 

+6.7 

1.1 

2.4 

2.8 

7 

700-845 

4.59 X 10' 

0.3 

0.5 

6 9 

1.3 

2.0 

2.8 

8 

845-1005 

1.46 X 10' 

0.6 

0.5 

+7.4 
-7 2 

1.4 

2.2 

2.8 

9 

1005-1195 

4.61 

0.9 

0.7 

+7.8 

+8:§ 

1.4 

2.9 

2.8 

10 

1195-1410 

1.39 

1.6 

0.7 

1.5 

3.6 

2.8 

11 

1410-1665 

4.14 X 10-' 

2.4 

0.6 

+^:i 

1.5 

3.6 

2.8 

12 

1665-1960 

1.11 X 10-' 

4.4 

0.6 


1.8 

3.4 

2.8 

13 

1960-5000 

4.33 X 10-3 

7.4 

0.7 

+i.9 

-9.3 

2.0 

3.0 

2.8 


Table 39; Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncer¬ 
tainties in the measurement. The events are selected using the inclusive analysis cuts and > 2, as well as 

Pj^ > 250 GeV. All other details are as for table 4. 
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Bin 

Bin edges [GeV] 

dcr/d(S/?“"‘™‘) [fb/GeV] 

« [%] 

[%] 

MJES [%] 

mjer [%] 

^unfold [%] 

^lumi 

1 

120-185 

1.38 X 10^ 

0.9 

2.1 

+7.1 

-7 1 

5.1 

8.3 

2.8 

2 

185-270 

3.96 X 10' 

0.5 

1.0 

+6.8 

7 0 

2.8 

5.7 

2.8 

3 

270-365 

5.95 X 10' 

0.4 

0.8 

+6.3 
-7 0 

1.8 

4.2 

2.8 

4 

365^65 

5.75 X 10' 

0.3 

0.8 

+6.6 

6 9 

1.6 

3.5 

2.8 

5 

465-575 

8.69 X 10' 

0.3 

0.5 

+6:4 

6 5 

1.3 

3.7 

2.8 

6 

575-700 

6.98 X 10' 

0.3 

0.5 

+6.9 

1.2 

3.7 

2.8 

7 

700-845 

3.85 X 10' 

0.3 

0.5 

-M 

-7 0 

1.3 

3.5 

2.8 

8 

845-1005 

1.45 X 10' 

0.6 

0.5 

+7.4 
-7 2 

1.3 

2.8 

2.8 

9 

1005-1195 

4.58 

0.9 

0.8 

+7.8 

1.4 

3.0 

2.8 

10 

1195-1410 

1.38 

1.6 

0.7 

+8.^ 

-7 8 

1.5 

3.6 

2.8 

11 

1410-1665 

4.11 X 10-' 

2.4 

0.6 

+8.8 

7 9 

1.6 

3.6 

2.8 

12 

1665-1960 

1.10 X 10-' 

4.4 

0.6 

+40 

1.8 

3.4 

2.8 

13 

1960-5000 

4.30 X 10-3 

7.5 

0.7 

+i.9 

-9.2 

2.0 

3.0 

2.8 

Table 40; Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncer¬ 
tainties in the measurement. The events are selected using the inclusive analysis cuts and > 2, as well as 

Pj^ > 400 GeV. All other details are as for table 4. 


Bin 

Bin edges [GeV] 

dcr/d(2/?“"‘™‘) [fb/GeV] 

« [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold [^] 

^lumi 

1 

120-185 

1.12 

3.4 

3.7 

+S1 

-7 5 

3.6 

20.8 

2.8 

2 

185-270 

2.62 

1.7 

1.9 

+7.8 

7 3 

3.3 

8.7 

2.8 

3 

270-365 

3.75 

1.3 

1.4 

+7.6 

2.5 

3.6 

2.8 

4 

365-465 

5.76 

1.1 

1.3 

8 7 

2.2 

3.5 

2.8 

5 

465-575 

7.27 

0.9 

1.0 

M 
-11 

2.2 

5.1 

2.8 

6 

575-700 

9.64 

0.8 

0.8 

+7.8 
-7 2 

2.1 

5.5 

2.8 

7 

700-845 

9.83 

0.6 

0.6 

+7.6 
-7 3 

1.7 

4.5 

2.8 

8 

845-1005 

6.31 

0.9 

0.7 

+7.6 
-7 4 

1.4 

3.7 

2.8 

9 

1005-1195 

3.99 

1.0 

0.7 

+7.6 

1.3 

3.5 

2.8 

10 

1195-1410 

1.37 

1.6 

0.7 

+i.^ 

11 

1.3 

3.9 

2.8 

11 

1410-1665 

4.07 X 10-' 

2.5 

0.6 

-7 9 

1.4 

3.7 

2.8 

12 

1665-1960 

1.09 X 10-' 

4.4 

0.6 

+8.8 

1.7 

3.4 

2.8 

13 

1960-5000 

4.26 X 10-^ 

7.5 

0.7 

+i.8 

-9.0 

1.9 

3.0 

2.8 


Table 41; Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncer¬ 

tainties in the measurement. The events are selected using the inclusive analysis cuts and Aj/™ > 2, as well as 

Pj^ > 550 GeV. All other details are as for table 4. 
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Bin 

Bin edges [GeV] 

dcr/d(2/?“"«“>) [fb/GeV] 

Cf [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi 

1 

120-190 

5.25 X 10 ^ 

0.7 

3.2 

10 1 

5.4 

8.0 

2.8 

2 

190-285 

5.85 X 10 ^ 

0.4 

1.3 

+ 9.7 
-9 7 

4.6 

6.5 

2.8 

3 

285-385 

1.77 X 10 ^ 

0.7 

1.1 

+ 8.8 
-8 0 

2.2 

4.5 

2.8 

4 

385^90 

5.1 X 10 ^ 

1.0 

0.9 

+ 8.0 

6 4 

1.2 

3.4 

2.8 

5 

490-605 

1.55 X 10 ^ 

1.1 

0.8 

+ 7.1 

6 2 

1.2 

2.3 

2.8 

6 

605-735 

4.7 X 10 ^ 

0.6 

0.9 

+7:t 

6 6 

1.2 

2.4 

2.8 

7 

735-880 

1.45 X 10 ' 

0.5 

0.8 

+ 7.6 
-7 3 

1.2 

3.5 

2.8 

8 

880-1040 

4.49 

0.9 

1.0 

+ 8.1 
-7 8 

1.2 

3.6 

2.8 

9 

1040-1225 

1.44 

1.4 

1.2 

+ 8.6 

-E 

8 9 

1.2 

4.0 

2.8 

10 

1225-1430 

4.69 X 10 ' 

2.9 

1.1 

1.2 

5.1 

2.8 

11 

1430-1655 

1.52 X 10 ^' 

4.5 

1.3 

1.2 

4.5 

2.8 

12 

1655-1905 

4.60 X 10-2 

7.4 

1.3 

1.5 

4.5 

2.8 

13 

1905-5000 

2.01 X 10-2 

10.7 

1.2 

+ 10 1 
- 10.7 

1.6 

6.8 

2.8 


Table 42: Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncertain¬ 
ties in the measurement. The events are selected using the inclusive analysis cuts and > 3. All other details 

are as for table 4. 


Bin 

Bin edges [ GeV ] i 

dcr / d ( 2 / 7 “'"'^) [ fb / GeV ] 

« [%] 

m 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi [^] 

1 

120-190 

1.86 X 10 ^ 

3.0 

2.4 

+ 9.1 

5 3 

2.0 

4.1 

2.8 

2 

190-285 

3.54 X 10 ^ 

1.3 

1.2 

+ 8.6 
-7 0 

2.0 

5.4 

2.8 

3 

285-385 

4.27 X 10 ^ 

1.2 

1.3 

+ 8.0 
-7 4 

2.0 

7.8 

2.8 

4 

385-490 

3.77 X 10 ^ 

0.9 

0.8 

+ 7.7 
-1 1 

1.9 

5.2 

2.8 

5 

490-605 

1.56 X 10 ^ 

1.1 

0.8 

+ 7 :i 

6 6 

1.5 

2.8 

2.8 

6 

605-735 

4.75 X 10 ' 

0.6 

0.9 

+ 7:2 

6 7 

1.4 

2.6 

2.8 

7 

735-880 

1.46 X 10 ' 

0.5 

0.8 

+ 7.6 
-7 4 

1.4 

3.5 

2.8 

8 

880-1040 

4.53 

0.9 

0.9 

+ 7.9 
-8 0 

1.4 

3.6 

2.8 

9 

1040-1225 

1.45 

1.4 

1.1 

+ 8.6 

1.4 

4.0 

2.8 

10 

1225-1430 

4.73 X 10 -' 

2.9 

1.1 

+i:i 

1.4 

5.1 

2.8 

11 

1430-1655 

1.53 X 10 -' 

4.5 

1.3 

+i;i 

1.5 

4.5 

2.8 

12 

1655-1905 

4.63 X 10-2 

7.4 

1.3 

+ 9.5 

1.9 

4.5 

2.8 

13 

1905-5000 

2.02 X 10-2 

10.6 

1.2 

1 + 1 

2.1 

6.8 

2.8 


Table 43: Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncer¬ 
tainties in the measurement. The events are selected using the inclusive analysis cuts and > 3, as well as 

Pj ^ > 250 GeV. All other details are as for table 4. 
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Bin 

Bin edges [GeV] i 

do-/d(S/7;;,™‘-'"‘) [fb/GeV] [%] 

[%] 

MJES [%] 

mjer [%] 

^unfold [^] 

^lumi 

1 

120-190 

3.81 

1.6 

4.0 

+9.3 

1.4 

3.6 

2.8 

2 

190-285 

1.35 X 10* 

0.8 

1.5 

+§.f 

1.4 

4.6 

2.8 

3 

285-385 

2.17 X 10' 

0.6 

1.1 

-1 3 

1.4 

6.1 

2.8 

4 

385^90 

2.12 X 10' 

0.6 

1.5 

+7.1 

6 9 

1.4 

7.9 

2.8 

5 

490-605 

3.56 X 10' 

0.4 

0.8 

+ 7:3 

6 9 

1.4 

5.1 

2.8 

6 

605-735 

2.46 X 10' 

0.5 

0.8 

+6.9 
-7 0 

1.4 

3.5 

2.8 

7 

735-880 

1.32 X 10' 

0.6 

0.9 

+ 1.2 

1.5 

3.0 

2.8 

8 

880-1040 

4.49 

0.9 

1.0 

-7 9 

1.9 

3.1 

2.8 

9 

1040-1225 

1.43 

1.4 

1.1 

+ 8.3 

2.0 

3.9 

2.8 

10 

1225-1430 

4.69 X 10-' 

2.9 

1.1 

+i.§ 

2.0 

5.1 

2.8 

11 

1430-1655 

1.51 X 10-' 

4.5 

1.3 

+i:l 

2.0 

4.5 

2.8 

12 

1655-1905 

4.59 X 10-2 

7.4 

1.3 

+^.4 

-9 1 

2.0 

4.5 

2.8 

13 

1905-5000 

2.00 X 10-2 

10.7 

1.2 

+ 9.9 

=ms 

1.7 

6.8 

2.8 

Table 44; Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncer¬ 
tainties in the measurement. The events are selected using the inclusive analysis cuts and > 3, as well as 


( 1 ) 


Pi 


> 400 GeV. All other details are as for table 4. 


Bin Bin edges [GeV] dcr/d(£j?^™‘”‘) [fb/GeV] [%] 


1 

120-190 

1.63 X 10-' 

7.7 

10.2 

2 

190-285 

4.76 X 10-' 

4.0 

4.0 

3 

285-385 

9.39 X 10 ' 

2.6 

2.8 

4 

385^90 

1.79 

2.2 

2.8 

5 

490-605 

2.21 

1.8 

2.1 

6 

605-735 

3.43 

1.2 

1.3 

7 

735-880 

2.86 

1.2 

1.0 

8 

880-1040 

2.10 

1.4 

1.3 

9 

1040-1225 

1.34 

1.4 

1.1 

10 

1225-1430 

4.64 X 10-' 

2.9 

1.1 

11 

1430-1655 

1.50 X 10-' 

4.5 

1.3 

12 

1655-1905 

4.54 X 10-2 

7.4 

1.3 

13 

1905-5000 

1.98 X 10-2 

10.7 

1.2 


mjes [%] 


+ 12.8 

-+^ 3^6 


-9.8 

+9.6 

-10^ 

+8.^ 

-9.1 

+8.3 

+ll 

- 8.0 

+ 8.1 


-8.9 

+9.7 

-10.5 


mjer [%] 
5.6 
5.6 

5.6 
5.5 
4.9 
3.0 
2.2 
2.1 
1.8 

1.7 
1.7 
1.7 
1.7 


llunfold [%] 
+8 
3.8 

4.2 

5.3 
6.2 

7.3 
6.1 
3.8 

3.7 
5.1 
4.6 
4.5 

6.8 


lllumi [^] 

Z8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 


Table 45; Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncer¬ 

tainties in the measurement. The events are selected using the inclusive analysis cuts and Aj/™ > 3, as well as 

^ > 550 GeV. All other details are as for table 4. 
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Bin 

Bin edges [GeV] < 

dcr/d(E/7“"“^) [fb/GeV] 

« [%] 

« [%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi [^] 

1 

120-190 

1.56 X 10^ 

1.1 

6.6 

+TU3 

-8 7 

4.2 

18.1 

2.8 

2 

190-285 

1.72 X 10^ 

0.8 

2.7 

+ 10.1 
-8 8 

4.2 

11.0 

2.8 

3 

285-385 

4.81 X 10^ 

1.3 

2.0 

+ 9.5 

3.7 

5.6 

2.8 

4 

385-490 

1.26 X 10^ 

2.1 

2.2 

+1.7 
-8 2 

2.4 

4.0 

2.8 

5 

490-605 

3.72 X 10' 

2.7 

1.7 

+ 6.9 
-7 9 

1.9 

3.2 

2.8 

6 

605-730 

1.09 X 10' 

1.8 

1.9 

+ 6.4 

1.8 

2.9 

2.8 

7 

730-865 

3.18 

1.1 

1.9 

+ 7.1 

1.8 

3.5 

2.8 

8 

865-1010 

1.07 

2.0 

2.6 

+U 

-9 8 

1.8 

5.5 

2.8 

9 

1010-1170 

3.49 X 10-' 

3.7 

3.2 

+ 8.9 
-12 0 

1.8 

5.9 

2.8 

10 

1170-1340 

1.29 X 10-' 

5.6 

3.6 

+ 9.1 

11 5 

1.8 

4.9 

2.8 

11 

1340-1525 

4.99 X 10-2 

9.8 

3.3 

+ 9.1 
-11 1 

1.8 

4.5 

2.8 

12 

1525-5000 

1.25 X 10-2 

12.0 

2.6 

+ 11.7 
-11.8 

1.8 

4.5 

2.8 


Table 46: Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncertain¬ 

ties in the measurement. The events are selected using the inclusive analysis cuts and Ai/™™ > 4. All other details 
are as for table 4. 


Bin 

Bin edges [GeV] i 

dcr/d(E/?“"«“') [fb/GeV] 

Cf [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi 

1 

120-190 

3.57 X 10' 

7.3 

5.7 

+3.6 

6 0 

2.0 

5.4 

2.8 

2 

190-285 

6.91 X 10' 

3.0 

2.8 

+7:3 

-7 7 

2.0 

6.1 

2.8 

3 

285-385 

9.82 X 10' 

3.1 

2.4 

+8.1 
-8 5 

2.0 

7.0 

2.8 

4 

385^90 

9.06 X 10' 

2.0 

1.7 

+7.7 

8 2 

2.0 

4.4 

2.8 

5 

490-605 

3.75 X 10' 

2.6 

1.8 

+r.i 

2.0 

3.1 

2.8 

6 

605-730 

1.1 X 10' 

1.7 

1.8 

+?:i 

2.0 

3.5 

2.8 

7 

730-865 

3.21 

1.1 

1.9 

+ 1.1 

-9 1 

2.0 

5.0 

2.8 

8 

865-1010 

1.09 

2.0 

2.5 

+ 8.5 
-10 2 

2.0 

5.8 

2.8 

9 

1010-1170 

3.53 X 10-' 

3.6 

3.2 

+9.2 

2.0 

6.5 

2.8 

10 

1170-1340 

1.30 X 10-' 

5.6 

3.5 

-11 0 

2.0 

6.3 

2.8 

11 

1340-1525 

5.03 X 10-2 

9.7 

3.3 

+^.5 

1 1 9 

2.0 

4.9 

2.8 

12 

1525-5000 

1.27 X 10-2 

11.9 

2.5 

+ 12.0 
- 19.1 

2.0 

4.1 

2.8 


Table 47: Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncer¬ 

tainties in the measurement. The events are selected using the inclusive analysis cuts and Aj/™’' > 4, as well as 

Pj^ > 250 GeV. All other details are as for table 4. 
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Bin 

Bin edges [GeV] < 

dcr/d(E/7“"‘™‘) [fb/GeV] 

KT [%] 

« [%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi [^] 

1 

120-190 

2.99 X 10-‘ 

5.3 

13.9 

+9.6 

6 3 

2.3 

5.6 

2.8 

2 

190-285 

1.53 

2.3 

4.9 

+g:? 

2.3 

5.6 

2.8 

3 

285-385 

2.72 

1.7 

3.5 

-9 9 

2.3 

6.3 

2.8 

4 

385-490 

2.92 

1.5 

4.3 

+9.5 

9 8 

2.3 

8.2 

2.8 

5 

490-605 

6.36 

1.0 

1.7 

+8.S 

-8 9 

2.3 

7.9 

2.8 

6 

605-730 

4.84 

1.1 

1.9 

+7.8 

2.3 

4.9 

2.8 

7 

730-865 

2.80 

1.2 

1.8 


2.3 

3.7 

2.8 

8 

865-1010 

1.08 

2.0 

2.6 

+i:? 

-9 7 

2.3 

4.2 

2.8 

9 

1010-1170 

3.49 X 10-' 

3.7 

3.3 

+9.6 
-11 4 

2.3 

6.0 

2.8 

10 

1170-1340 

1.29 X 10-' 

5.6 

3.5 

+9.6 
-11 9 

2.3 

6.2 

2.8 

11 

1340-1525 

4.99 X 10-2 

9.8 

3.3 

+^.4 
-11 7 

2.3 

4.9 

2.8 

12 

1525-5000 

1.25 X 10-2 

12.0 

2.5 

+ 11.1 
-P..7 

2.3 

4.1 

2.8 


Table 48: Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncer¬ 

tainties in the measurement. The events are selected using the inclusive analysis cuts and Aj/™’' > 4, as well as 

Pj^ > 400 GeV. All other details are as for table 4. 


Bin 

Bin edges [GeV] 

dcr/djEp^"'™') [fb/GeV] 

« [%] 

[%] 

mjes [%] 

mjer [%] 

^unfold [%] 

^lumi [^] 

1 

120-190 

4.02 X 

10-* 

103.7 

141.7 

+1X0 

1 8 6 

3.2 

4.8 

2.8 

2 

190-285 

1.92 X 

10-2 

20.4 

30.3 

3.2 

17.8 

2.8 

3 

285-385 

4.80 X 

10-2 

11.4 

14.6 

+ 16.1 
-18 2 

3.2 

23.4 

2.8 

4 

385-490 

1.51 X 

10-' 

6.8 

9.1 

+ 16.2 
-17 0 

3.2 

23.9 

2.8 

5 

490-605 

1.42 X 

10-' 

6.6 

8.0 

+ 15.2 

15 6 

3.2 

21.6 

2.8 

6 

605-730 

3.53 X 

10-' 

3.6 

3.7 

+ 12.5 
-124 

3.2 

14.8 

2.8 

7 

730-865 

3.55 X 

10-' 

3.8 

2.7 

+ 10.2 
Ti.V 

-11 7 

3.2 

10.7 

2.8 

8 

865-1010 

3.75 X 

10-' 

2.9 

3.4 

3.2 

6.6 

2.8 

9 

1010-1170 

3.03 X 

10-' 

4.0 

3.1 

3.2 

5.0 

2.8 

10 

1170-1340 

1.28 X 

10-' 

5.6 

3.5 

3.2 

4.9 

2.8 

11 

1340-1525 

4.95 X 

10-2 

9.8 

3.3 

+10.6 
-12 0 

3.2 

4.9 

2.8 

12 

1525-5000 

1.24 X 

10-2 

12.1 

2.6 

+ 10.2 

-12.6 

3.2 

3.7 

2.8 


Table 49: Measured differential four-jet cross section for R - 0.4 jets, in bins of along with the uncer¬ 

tainties in the measurement. The events are selected using the inclusive analysis cuts and Aj/™’' > 4, as well as 

Pj^ > 550 GeV. All other details are as for table 4. 
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